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ABSTRACT 

We use 317,000 emission-line galaxies from the Sloan Digital Sky Survey to investigate line-ratio se¬ 
lection of active galactic nuclei (AGNs). In particular, we demonstrate that “star formation dilution” 
by Hii regions causes a significant bias against AGN selection in low-mass, blue, star-forming, disk- 
dominated galaxies. This bias is responsible for the observed preference of AGNs among high-mass, 
green, moderately star-forming, bulge-dominated hosts. We account for the bias and simulate the in¬ 
trinsic population of emission-line AGNs using a physically-motivated Eddington ratio distribution, in¬ 
trinsic AGN narrow line region line ratios, a luminosity-dependent Lboi/A[Oiii] bolometric correction, 
and the observed Mbh — <7 relation. These simulations indicate that, in massive (log(M*/MQ) > 10) 
galaxies, AGN accretion is correlated with specific star formation rate but is otherwise uniform with 
stellar mass. There is some hint of lower black hole occupation in low-mass (log(M*/M 0 ) < 10) 
hosts, although our modeling is limited by uncertainties in measuring and interpreting the velocity 
dispersions of low-mass galaxies. The presence of star formation dilution means that AGNs contribute 
little to the observed strong optical emission lines (e.g., [Oiil] and Ha) in low-mass and star-forming 
hosts. However the AGN population recovered by our modeling indicates that feedback by typical 
(low- to moderate-accretion) low-redshift AGNs has nearly uniform efficiency at all stellar masses, star 
formation rates, and morphologies. Taken together, our characterization of the observational bias and 
resultant AGN occupation function suggest that AGNs are unlikely to be the dominant source of star 
formation quenching in galaxies, but instead are fueled by the same gas which drives star formation 
activity. 

Subject headings: galaxies: active - galaxies: nuclei - galaxies: Seyfert - quasars: emission lines - 
galaxies: evolution 
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1. INTRODUCTION 

The observed correlations between the mass of 
a galaxy’s bulge and the m ass of its supermassive 
black hole (S M BH) fe. g.. iMaeorrian et al.l 1R98; 

iFerrare se &: M errittI 2000 : 

iKormendv fc Hoi 120131 1 im¬ 
ply that galaxy growth via star formation (SF) must 
have a corresponding period of SMBH growth in the 
active galactic nucleus (AGN) phase. Yet the details 
that couple AGN-galaxy coevolution remain myste- 
rious. Theoretical f rameworks invoke mergers (e.g., 
Sanders et al.l I1 988I |Di Matteo. S nringel. fc Hernouistl 
20051 : iHopkins et al.l l2006l 20(^ . “ violent disk 

instabilities (jPekel et al.l 120091 : iBournaud et al.l 
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1201 H iGabor fc BournaudI 1201311. or secular pro¬ 
cesses in gas-rich d i sks ( e. g.. iShlosman et al.l [T989I 
iHopkins fc Hernaiiistl 120091 iHookins et al.l 201411 to 
simul taneously (or near-simultaneously, e.g.. iWild et al.l 
1201011 fuel rapid star formation and drive gas inward 
to power an AGN. Many simulations also invoke AGN 
“feedback” which quenches s tar formation through gas 


blowout bv radiative winds (ISilk & Rees 

119981: lEabiaul 

120021 iDi Matteo. Sorb 

igel. & Hernouistl 

200511 or gas 

heating by radio jets ( 

Groton et al.1 1200fi 

). Both AGN 


feedback and coupled AGN-SF fueling leave an imprint 
on the observed properties (e.g., color, star formation 
rate, and morphology) of AGN host galaxies. This sce¬ 
nario means the physical processes behind AGN-galaxy 
coevolution can be revealed by the simple question: Are 

AGN host galaxies spe cial? _ 

In seminal work iKauffmann et all (|2003fl used 
^120,0 00 galaxies from the Sloan Digital Sky Survey 
ISOSS. lYork et alll2000ll to demonstrate that AGNs are 
most frequently observed in massive bulge-dominated 
galaxies, with more luminous AGNs fo und in host galax¬ 
ies wi th more recent star formation. IKauffmann et ahl 
(|2003[1 . along with subseque n t wor k reaching similar 
concl usions (iHeckman et al.l 12004 iSchawinski et ahl 
[Moll . selected AGNs using the characteristic optical 
signature of the AGN narrow line region (NLR): AGN 
emission produces higher ratios of partially ionized 
forbidd en lines comp ar ed to Balmer recombination 
lines (jBaldwin et al'l 119811 : Veilleux fc OsterbrockI 
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pies 

2007 

witn A-rav iiiNana 
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20081 iSilverman et al. 

2008 

iSchawinski et al.l 

20091 iHickox et all 120091 


Kocevski et all 120091 : iHaggard et al.l l20i nli. in frared 
(jHickox et al.ll2009[ i. and radio (jBest et al.l 1200,1) AGN 


selection found a similar observed preference for massive 
and bulge-dominated hosts, with infrared X-ray 
radio corresponding to a sequence in d eclining star 
forma tion rate (SFR) in the host galaxy (|Hickox et alJ 
I2009H . Taken together these observations implied a 
critical stellar mass for accretion onto an AGN, with 
indirect evidence for AGN feedback in the green/red 
colors of X-ray/radio AGN host galaxies. 

More recent studies, however, have suggested that 
the apparent green-valley prefere nce of X-ray AGNs is 
purely a relic o f selection effects (|Silverman et al.l [20091 
IXue et ;ini201fl[ l . There is als o a mass bias for AGN selec¬ 
tion at a fixed accretion rate (|Aird et al.ll 2012 ll . most eas¬ 
ily understand with the unitless Eddington ratio, AEdd = 
Aboi/EEdd- The bolometric luminosity is driven by mass 
accretion, Lboi = rjMc^, and the Eddington luminosity 
is defined as AEdd = (1-3 x 10^® erg s“^)MBH/Af 0 , so 
that Aboi/AEdd — 0.44?7(M/Gyr“^)/MBH- The radiative 
efficiency 77 ^ 0.1 at m = M/[L^dd/iv^'^)] ^ 10 “^, while 
lower accretion rates lead to radiatively i nefficient ac¬ 
cretio n and log( 77 ) ^ —1 -I- 0.5 log(m/10“^) (|Xie fc YuanI 
I 2 OI 2 D . The Eddington ratio better describes the mode of 
AGN accretion and feedback than does luminosity, as the 
AGN accretion flow changes from a geometrically thin 
accretion disk with radiative winds at AEdd ^ 0.01 to a 
radiatively inefficient accretion flow with stronger radio 
jets at Avdd 5. 0.01 (iNaravan fc McGlinto^ 120081 iHol 
120081 iTrurnp et al.l l2011al IkRckman fc BestI 120141 1. Be¬ 
cause Mbh is correlated with galaxy bulge mass Mbuige, 
which is in turn correlated with total stellar mass M*, a 
given AGN flux limit means it is easier to find lower-AEdd 
AGNs in high-mass galaxies than in low-mass galaxies. 
The AEdd function for AGNs is likely to steeply decline 
with increasing accretion rate, since there are many more 
weakly-accreting t han rapid ly-accreting AGNs in the lo¬ 
cal universe (e.g., iHol I 2 OO 8 D . In a flux-limited sample 
this results in a steep bias towards finding more AGNs 
i n massive hosts. 

lAird et HI (j2012( l reported that the X-ray AGN host 
occupation fraction is consistent with a uniform AGN 
Eddington ratio distribution across all host galaxy stel¬ 
lar masses. Thus, after controlling for selection bias, 
AGN host galaxies are not special in terms of stellar 
mass. It remains unclear if AGN host galaxies are spe¬ 
cial in SFR or morphology. Recent work combining far- 
infrared observations with X-ray selection demonstrates 
that AGNs are most common in star-forming galaxies 


(iMullanev et al.ll2012l:lHarrison et al.ll2012 

Juneau et alJ 

120131 IRosario et al.l I2013allbl: IGhen et al. 

2 ?nib iris 


result is generally similar to the connection between 
AGN luminosity and host gal axy SFR seen in the ini - 
tial line-ratio AG1N0 sample of iKauffmann et al.l (120031 1 
and in subsequent studies of broad-line AGN hosts 
(|TrumD et al.ll201^ iMatsuoka et al.l[2014ll . In terms of 

® Throughout the paper we refer to AGNs selected by Seyfert- 
like optical line ratios as “line-ratio AGNs.” 


host morphology, most recent observations demonstrate 
that AGNs are not p r eferentially fueled by ma jor mergers 


(IGrogin et al.l 120051 

Gabor et al.l 120091 ICisternas et al.l 

120111 iKocevski et al. 

120121). althou 

gh mergers may play 
it are verv luminous 

a role in fueling nearby AGNs thi 

and/or obscured (IKoss et al.l 12010 

: lEllison et al.l 120111 

iTreister et al.l 120121 iTrumol l2013b 
mixed on whether violent disk ins 
nificant role in AGN fueling, wit 
(|Bournaud et al.ll2012n and against 
So while X-ray AGN hosts do not 
in their stellar mass, it remains unc 
special in SFR or morphology. And 
the larger population of line-ratio 
N-ray picture. 

In this work we use line-ratio AG 

). Observations are 
tabilities play a sig- 
h evidence both for 
(|TrumD et alJl2014ll. 

appear to be special 
ler debate if they are 
it is still unclear how 
AGNs connect to this 

M selection to recover 


the intrinsic AGN occupation fraction with host galaxy 
properties. For this purpose we use a very large sample 
of >300,0 00 galaxies from the SPS S, described in Section 
2. Just as IKauffmann et all (1200311 originally found, Sec¬ 
tion 3 demonstrates that line-ratio AGNs are observed to 
be most common in massive, concentrated galaxies with 
intermediate colors and specific star formation rates (sS- 
FRs). The selection bias for line-ratio AGNs is a func¬ 
tion of the contrast between AGN and galaxy emission 
lines rather than a simple flux limit, and in Section 4 
we quantify this bias with host galaxy properties. In 
Section 5 we model the intrinsic AGN population as a 
function of galaxy properties, beginning with a uniform 
Eddington ratio distribution, then allowing Eddington 
ratio and black hole occupation to vary with galaxy prop¬ 
erties. Section 6 describes the implications of the mod¬ 
eled AGN occupation function for galaxy emission-line 
measurements, the connection between AGN accretion 
and galaxy SFR, AGN feedback, and black hole seed for¬ 
mation. Throughout the paper we use a basic AGDM 
cosmology with h = 0.70, (Im = 0.3, Ha = 0.7. 

2. OBSERVATIONAL DATA 

Our galaxy sample i s drawn from th e Sloan Digi¬ 
tal Sky Survey (SDSS, lYork et ^ I200C ). SPSS data 
inclu de ugriz broad-band photometry ( Fukugita et al.l 
Il990fl and spectroscopy with R ~ 2000 over 3800 < 
A < 9200A (ISmee et ahllMll . taken using a 2.5-m tele¬ 
scope at Apache Point Observatory (iGunn et al.ll2006ll . 
The galaxy sampl e in this work is drawn from the SDSS 
Data Release 10 (|Ahn et al.l 1201411 . which covers a total 
sky area of 14,555 deg^. We used the public SkyServer 
SQL servei0 to select a parent sample of 317,192 galax¬ 
ies (class=GALAXY) with r < 17.77, spectroscopic red- 
shifts in the range 0.01 < z < 0.1, and line flux errors of 
0 < (Tiine < 10“^® erg s“^ cm“^ (to remove bad spectra). 
The r-magnitude limit represents t he main spectroscop y 
limit of the SDSS galaxy survey (iStrauss et al.l [200^ . 
The limits in the line flux errors effectively remove ob¬ 
jects with problematic spectra (artifacts and bad pixels) 
in the line regions. These spectral issues are uncorre¬ 
lated with galaxy properties and are essentially random 
occurrences. 

The redshift range of 0.01 < z < 0.1 is designed 
to include a broad range of galaxy stellar masses with¬ 
out significant evolution in galaxy properties. The up- 

® http://skyserver.sdss3.org/drl0/en/tools/search/sql.aspx 
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per limit in redshift corresponds roughly to a mass 
limit of log(M*/M©) > 10.7 with the SDSS magni¬ 
tude limit, and the lower redshift limit adds lower- 
mass galaxies to log(M*/M 0 ) > 8 . In Section 4.1 we 
demonstrate that the fractional coverage of the SDSS 
fiber varies by < 20 % over the stellar mass range of 
the full sample. There is also minimal galaxy evolu¬ 
tion within this redshi ft range: for example, the cosmic 
SFR evolves by <25% llWyder et al.ll2005tlCucciati et al.l 
120121: iMadau fc Dickinsonil2014l l. Thus we do not expect 
a significant change in AGN fraction with redshift driven 
by changing galaxy properties. 

SDSS galaxies have redshifts and line measurements 
computed using the idlspec2d sof tware: details of thi s 
automated processing are given by iBolton et al.l (I2012D . 
In brief, galaxy classifications and redshifts are computed 
by finding the best-fit (minimum-y^) template, with tem¬ 
plates constructed from sets of eigenspectra derived from 
a principal component analysis of training spectra. Emis¬ 
sion line fluxes are computed using the best-fit velocity 
dispersion template as the continuum, which implicitly 
corrects for stellar absorption. Each line is fit as a Gaus¬ 
sian, with line width fixed across the Balmer lines, and a 
separate fixed line width for the group of all other emis¬ 
sion lines visible in our redshift range. 

Stellar velocity dispersion is computed by fit- 
ting a set of stellar tem plates from the library of 
iPrugniel fc SonbiranI (j2001h . degrading the templates to 
match the instrumental broadening (~70 km s“^) and 
minimizing over a grid of 0 < ct < 850 km s“^ with 
25 km s“^ bins. Since these velocity dispersions are fit 
using the stellar absorption lines, they are unaffected by 
AGN emission. Low-mass galaxies (log(M*) < lO^^M©) 
tend to have measured dispersions less than the instru¬ 
mental resolution: we include dispersions for these galax - 
ies using survival analysis fe.g.. lFeigelson fc NelsoDlll985ll 
on their l-tr upper limits. The lack of reliable velocity 
dispersions in low-mass galaxies is an issue we return to 
in Sections 5 and 6. 

We create sam ples of galaxies with well-measured line 
ratios, following I.Tunean et al.l (j2014l l and requiring the 
signal-to-noise ratio (S/N) of the line ratios to have 
(S'/A^) ratio > 3/-\/2. This is the equivalent of each line 
detected at {S/N)iine > 3, but has the advantage of 
including well-constrained line ratios where one line is 
bright and the other faint (for example, H/3 is frequently 
bright and [Oiii] weak in massive galaxies, and Ha is 
bright and [Nil] weak in low-metallicity galaxies). From 
the parent sample of 317,192 galaxies with r < 17.77 
and 0.01 < z < 0.1, 223,448 galaxies possess each of 
the [Oiil]A5007/H/3 and [Nii]A6584/Ha ratios measured 
at (5'/iy)ratin > 3/ \/2: these lines form the “BPT” 
([Baldwin et al.lHoMIl diagram and so we designate this 
galaxy set as the “BPT well-measured sample.” Sim¬ 
ilarly the [Oiii]A5007/H/3 and f[Sii]A6717-H 67311/Ha 
ratios define the “V087” (|Veilleux fc OsterbrockI [l9^ 
diagram, and the ‘W087 well-measured sample” in¬ 
cludes the 215,242 galaxies with these ratios measured 
at (S'/A')ratio > 3/-\/2. For most galaxies the (5'/A^)ratio 
threshold corresponds to a line flux limit of / > 1 x 
10“^® erg s“^ cm“^. Our line-ratio classifications gen¬ 
erally use only these well-measured samples, although 
we also discuss the likely AGN fraction among the full 


parent sample of galaxies (including those with poorly- 
measured line ratios) in Section 3.3. 

Unless otherwise specified, the emission-line luminosi¬ 
ties in each galaxy are du st-corrected using the measured 
Balmer decreme nt and a iCardell i et al.l ()1989H attenua¬ 
tion curve with lO’Donnelll ( 19941 ) coefficients. The ex¬ 
tinction corrections use an intrinsic (dust-free) Balmer 
decrement of Ha/H,5 = 2.86 for SF-dominated emission¬ 
lines (i.e., classified as Hii-dominated on the BPT or 
V087 diagrams) and Ha/H/3 = 3.1 for AGN emission 
(appropriate for narrow line region gas conditions, e.g., 
lOsterbrock fc Ferlandl 1200611 . Galaxies with observed 
Balmer decrements less than the assumed intrinsic decre¬ 
ments were assigned zero extinction. 

We characterize basic galaxy properties using stellar 
mass, color, specific star formation rate (sSFR), and con¬ 
centration. Stellar ma sses and specific star formatio n 
rates are estimated by iMontero-Dorta et al.l din oreo.l l , 
fitting the broad-band ugriz SDSS photometry with a 
grid of templates fr om the flexible stella r population syn¬ 
thesis co de (FSPS , iGonrov et al.l 1200911 . The templates 
assume a iKroupgil (|2001[ 1 initial mass function, allowing 
a wide range of galaxy formation times 2 < tage/Gyr < 
Tu{z), where Tu is the age of the universe at redshift z. 
The stellar mass and sSFR determination uses a broad¬ 
band fit for the dust attenuation (separate from the 
Balmer decrem ents used for du s t-corr ecti ng the emission 
lines), following iCharlot fc Falll (|200n[ l and lGalzetti et al.l 
(j200(tl by assuming ri for birth clouds and T 2 for old 
stars, with ti = 3r2, and fitting a grid of values span¬ 
ning 0 < r 2 < 0.75. The use of broad-band photometry 
to estimate stellar mass and sSFR minimizes potential 
AGN contamination. Our study uses only narrow-line 
AGNs (not quasars), which generally have weak contin¬ 
uum emission. The AGN emission lines also do not sig¬ 
nificantly affect the broad-band photometry, since almost 
all AGNs in the sample have emission line equivalent 
widths <100A. 

All galaxy colors are k-cor rected to z = 0.05 usiii g 
the kcorrect IDL software (|Blanton fc Row^ l2007fl . 
We characterize galaxy morphology using a concentra¬ 
tion parameter Cr dehned as the ratio of the radii 
containing 90% and 50% of the galaxy’s r-band light, 
Cr = Rein r /R h n r ■ Details on the m easurement of 
these iPetrosianI (|1976fl radii are given by iBlanton et al.l 
(|2001h . Galaxies with high concentration (Rgo.r/^^so.r ^ 
A 5) tend to be bulge-dominated, while galaxies with 
i? 90 ,r-/i? 5 o,r ^ 2.5 are typically disks. 

It is important to note that the sample of galaxies with 
well-detected emission lines occupies a slightly different 
parameter space in galaxy properties compared to the full 
SDSS parent sample. In Figure[l]we compare the color- 
mass, sSFR-mass, and concentration-mass diagrams for 
the full parent sample with the well-measured BPT and 
V087 samples. The restrictions on BPT and V087 line 
ratios result in nearly identical well-measured samples, 
with fewer red, concentrated, and low-sSFR galaxies. 
Since the galaxies removed from the well-measured sam¬ 
ple typically have S/Nune < 3, they tend to have weaker 
emission lines and so are significantly less likely to have 
AGNs. In Section 3.3 we devise a strategy that accounts 
for the small AGN fraction likely present in galaxies with 
poorly-measured line ratios. 
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All Galaxies BPT Ratias S/N>2.12 


BPT S/N>2.12 V087 S/N>2.12 



8 9 10 11 128 9 10 11 12 89 10 11 128 9 10 11 12 

Stellar Mass [log(M0)] Stellar Mass [log(MQ)] Stellar Mass [log(Mg)] Stellar Mass [log(Mg)] 


Fig. 1.— The distribution of galaxy properties for the full parent sample (left column, 317,192 galaxies) with the well-measured sample 
of galaxies having both BPT line ratios measured at (S/fV)ratio > 3/\/2 (second column, 223,448 galaxies). The two columns at right also 
display the fraction of the full sample with well-measured ((S/A'^)ratio > 3/\/2) BPT or V087 line ratios. The emission-line ratio S/N 
requirements result in fewer red, low-sSFR, concentrated galaxies. 


3. OBSERVED HOST PROPERTIES OF LINE-RATIO AGNS 

We begin by examining the observed AGN fractions 
among host galaxies of different mass, color, SFR, and 
morphology. Section 3.1 introduces our adopted meth¬ 
ods of line-ratio selection, including a method to sepa¬ 
rate AGNs and low-ionization narrow emission regions 
(LINERs). We use the well-measured sample in Section 
3.2 to determine the fraction of BPT and V087 AGNs 
with galaxy color, specific star formation rate, morphol¬ 
ogy, and stellar mass. In Section 3.3 we additionally 
provide a correction for the (small) fraction of AGNs 
missed among galaxies with poorly-measured line ratios, 
concluding with an estimate of the overall AGN fraction 
among the entire parent sample. 

3.1. Line-Ratio AGN Selection 

The most widely used metho ds for line-ratio class ifica- 
tion of A GNs use the “BPT” dBaldwin et al.lll98il) and 
“V087” (|Veillemc fc Osterbrocklll987[ ) diagrams. Each 
compares the ratios of forbidden lines of partially ion¬ 
ized metals with Balmer recombination lines of hydro¬ 
gen: [Oiii]A5007/H/3 vs [Nil]A6584/Ha for the BPT, and 
[Oiii]A5007/H^ vs [Sii]A(6717 + 6731)/Hq; for the V087 


diagram. These lines are also typically among the bright¬ 
est features in rest-frame optical spectra of galaxies, and 
the close wavelength range of each pair means that the 
ratios are largely insensitive to dust (assuming that for¬ 
bidden metal lines and hydrogen recombination lines are 
emitted in the same galaxy regions). Note that BPT 
and V087 classification is appropriate for use on narrow 
(i.e. <1000 km/s) emission lines only. Our sample does 
not include broad-line AGNs, which also exhibit broad 
(>1000 km/s) Ho: and H/3 emission lines in their spectra. 

Eigure [5] presents our well-measured galaxy samples 
in the BPT and V087 diagrams. Each galaxy in this 
sample has both of the relevant line ratios measured at 
(S'/A^)ratio > 3/-\/2 (at left, [Oiii]/H/3 and [Nil]/Ha, and 
at right, [Oiil]/H/3 and [SiiJ/Ha). Also shown in Figure 
[2] are several demarcation lines for classifying galaxies: 

1. BPT [l^uffmann et al.l (I2003D empirical AGN/SF 
division line (long-dashed red), 

log([Oiii]/H/3) = 1.3 -b 0.61/[log([Nii]/Ha) - 0.05] (1) 

2. V087 empirical AGN/SF division line defined for 
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BPT V087 



Fig. 2.— The BPT and V087 diagrams for line-ratio classifica¬ 
tion of AGNs and star-forming galaxies. Median line ratio errors 
are shown by the error bars in the lower left of each panel. We 
adopt the long-dashed red lines for AGN /SF classification : in th e 
BPT diagram this is the empirical line of IKauffmann et al. mnni), 
and in the V087 diagram this is an empirical lin e of our own 
constr uction. Also shown in each diagram are the IKewlev et al.l 
||200!I'I maximal starburst lines (short-dashed blue), and the green 
dotted line in the V Q87 diagram is the AGN/LINER division of 
IKewlev et al.l l|2Q06lh 

this work (long-dashed red), 

log([Oiii]/H/3) = 1.3 -h 0.48/[log([Sii]/Ha) - 0.10] (2) 

3. BPT IKewlev et al.l (l2001f) maximal starburst line 
(dashed blue), 

log([Oiii]/H/3) = 1.19-b0.61/[log([Nii]/Ha)-0.47] (3) 

4. VO87 IKewlev et'all (1200111 maximal starburst line 
(dashed blue), 

log([Oiii]/H/3) = 1.3 -b 0.72/[log([Sii]/Ha) - 0.32] (4) 

5. V087 lKewlev et ^ (|2006ll division between AGNs 
and LINERs (dotted green), 

log([Oiii]/H/3) = 0.76-b 1.89log([Sii]/Ha) (5) 

We adopt Equations 1 and 2 as our lines for AGN/SF 
classification in this work, along with Equation 5 for 
AGN/LINER separation. 

Accurate AGN selection lines are subject to consid¬ 
erable uncertainty, largely d ue to differen c es in inter¬ 
pretation. For example, the IKewlev et al.l (120011) lines 
(Equations 3 and 4) are the maximal line ratios for 
starbursts in photoionization models, and so can be 
treated as minimum AG N li nes. Galaxies betwee n the 
iKauffmann et al.l (|2003D and IKewlev et al.l ()2001l) lines 
(Equations 1 and 3) in the BPT diagram are often classi¬ 
fied as “composite galaxies” with emission-line contribu¬ 
tion from both SF and AGN. In general, com posites have 
highe r SFRs than AGN-dominated galaxies (ISalim et al.l 
I2007D . and the t wo categories exhib i t little difference in 
X-ray properties (lJuneau et al.ll20lil : iTrouille et al.ll201ll : 

iTrumrj et al.l l20Ilb[) . Therefore we assume that com¬ 
posite galaxies host similar AGNs but with more dilu¬ 
tion from Hii regions compare to galaxies with AGN- 
dominated line ratios. For this reason we include both 
composite and AGN-do minated galaxies lying above the 
IKauffmann et al.l ()2003D empirical AGN/SF line (Equa¬ 
tion I) in our initial BPT AGN selection. 



Fig. 3. — Distributions of dust-corrected L[Oiii] (left) and 
L[Oiii]/M* (right) for AGNs and LINERs classified on the V087 
diagram. The dotted lines show all objects, while the solid lines 
show only “well-separated” AGNs and LINERs with line ratios at 
least Icr from the AGN/LINER classification line (Equation 5). 
Vertical lines show empirical divisions between the two categories, 
with percentages of outliers indicated for the well-separated AGN 
(blue) and LINER (red) samples. The L[Oiii]/M* ratio, a proxy 
for Eddington ratio, better distinguishes the two categories than 
[Oiii] luminosity alone, and we adopt this ratio to separate AGNs 
and LINERs identified by the BPT diagram. 

We create a new empirical line in the V087 diagram 
(Equation 2) which is designed to be parallel to the 
star-forming galaxy seq uence in a similar fashion to the 
IKauffmann et all ()2003l) empirical BPT li ne. In particu¬ 
lar, ou r new V087 line contrasts with the IKewlev et al.] 
()2001l) maximal starburst line in the treatment of low- 
metallicity galaxies (the upper left in each panel of Figure 
HD. Equation 4 selects several AGNs from the symmetric 
tail of the low-metallicity star-forming galaxy sequence, 
while our Equation 2 avoids classifying these galaxies as 
AGNs. 

In the V087 diagram, the IKewlev et al.l (|2006f ) 
line (Equation 5) additionally separates AGNs from 
low-ionization narrow emission-line regions (LINERs, 
iHeckmanl 119801) . X-ray and radio observations sug gest 
that LINER galaxies are likely to host nuclear AGNs (iHol 
I2008L and references therein). However LINERs have 
much low er accretion rates than stronger-lined AGNs 
()Hcil I2009D and much of their line emission is proba¬ 
bly powered by e xtended ionization so urces rather than 
nuclear activity (lYan fc BlantonI I2012D . For these rea¬ 
sons we avoid LINERs in our sample and define “V087 
AGNs” as lying abo ve the new empirical AGN/SF line 
(Equation 2) and the IKewlev et al.l (l2006f ) AGN/LINER 
line (Equation 5) in the V087 diagram. 

Unlike the “V087 AGN” category, the BPT diagram 
does not distinguish between AGNs and LINERs. In¬ 
deed, the BPT diagram does not lend itself to a simi - 
lar LINER/AGN separation as the IKewlev et al.] (120061 ) 
line in the V087 diagram: while the [Sii] critical den¬ 
sity relates to the LINER/AGN transition, the [Nil] 
critic al density is a factor of ^ 30 higher and does 
not. ([Veilleux fc OsterbrockI (119871) originally suggested 
a simple AGN/LINER separation of [Oiii]/H/3 > 3, but 
AGNs and LINERs form a continu ous rather than bi- 
modal distribution in [Oiii]/H/3, e.g. iHo et al.ll2003l ) 

AGNs typically have stronger emission lines than LIN¬ 
ERs, and so we investigate the use of [Oiii] luminosity 
and the ratio L[Oiii]/M* to separate the two classes in 
Figure [3] The L[Oiii]/M* ratio is related to the Edding- 
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ton ratio Tboi/-^^Edd: for a bolometric correction fc[oiii] = 
-f^boi/i[Oiil] (the inverse of the fraction of total AGN 
light reprocessed as [O ill] emission) and a mass ratio 7 = 

I/[Oiii]/M* = i.sxio^s'mbh ’ 

ables Ajoiii] and 7 are not actually constants: more accu¬ 
rate relationships between L[Oiii]/M* and Eddington ra¬ 
tio are discussed in Section 5. However log(L[Oiil]/M*) 
does monotonically increase with log(Lboi/EEdd) with a 
nearly linear slope of ^0.5, making it a useful simple 
proxy.) 

Figure [3] compares L[Oiil] and L[Oiil]/M* of “well- 
separated” AGNs and LINERs with measured line ratios 
at least la from the AGN/LINER classification line. The 
left panel shows a significant tail of low-luminosity V087 
AGNs overlapping with the LINERs. This overlap is re¬ 
duced when the two populations are instead compared 
in L[Oiil]/M*. In the right panel of Figure[3l AGNs and 
LINERs with line ratios well-separated from the classifi¬ 
cation line are divided by: 

log(L[Oiil]/M*) > 29.5 [log(erg s'^ M/)^)] ( 6 ) 

Due to both the physical motivation and the good em¬ 
pirical separation, we use Equation 6 to define a “BPT 
AGN” category which is mostly free of LINERs. 

In summary, we identify two categories of AGNs using 
emission-line properties: 

• BPT AGN: Above the lKauffmann et al.l (|2003f) line 
log([Oiii]/H/3) > 1.3-H0.61/[log([Nii]/Ha)-0.05], 
and with log(L[Oiii]/M*) > 29.5 (units of erg s“^ 
and Mq): Equations 1 and 6. 

• V0 87 AGN: Abov e our em pirical AGN/SF line and 
the iKewlev et al.l (|2006f l AGN/LINER division, 
log([Oiii]/H/3) > 1.3 -h 0.48/[log([Sii]/Ha) - 0.10] 
and log([Oiii]/H/3) > 0.76 -I- 1.891og([Sii]/HQ!): 
Equations 2 and 5. 

Among each well-measured sample, 6.3% 
(14,175/223,448) of galaxies are classified as BPT 
AGNs and 4.3% (9210/215,242) are classified as V087 
AGNs. In the next section we demonstrate that both 
methods identify AGNs in similar kinds of galaxies, and 
the lower AGN fraction identified by the V087 is due to 
slightly lower overall selection efficiency. 

3.2. Observed AGN Fractions in Well-Measured 
Galaxies 

Figure [4] shows the observed BPT and V087 AGN 
fractions in bins of galaxy color, specific star forma¬ 
tion rate, and concentration versus stellar mass among 
galaxies of the well-measured sample. The AGN frac¬ 
tion is defined as the fraction of galaxies (with line ratio 
(S'/A^)ratio > 3/-\/2) in each bin that meet the criteria of 
the BPT AGN or V087 AGN categories defined in Sec¬ 
tion 3.1, and bins with less than 50 galaxies are left blank. 
In both the BPT and V087 diagrams, the AGN fraction 
is observed to be highest in massive (M* ^ lO^^'^M©), 
green, moderate-sSFR, bulge-dominated galaxies with 
Rbo.t/R 50 ,r > 2.5. There is also a small population of 
V087 AGNs in low-mass hosts which are not identified 
by the BPT AGN classification. We demonstrate in Sec¬ 
tion 5 that this result may be due to a bias of the BPT 


against AGNs with low-metallicity NLRs which remain 
detected in the V087 diagram. 

The maximal AGN fraction in massive green, 
moderate-sSFR, bulge-dominated galaxies led to previ¬ 
ous conclusions that AGN require massive host galax¬ 
ies a nd cause feedback that quenches star formatio n 
le.g.. iKauffmann et al.l 120031 iSchawinski et al.l 12010(1 . 
iXue et al.l (l2010f) demonstrated that the similar massive 
green-valley host preference observed for X-ray AGNs is 
caused solely by selection effects, and lAird et id] (|2012f l 
additionally showed that the intrinsic AGN distribution 
is actually uniform at all stellar masses. In Sections 4 and 
5 we similarly investigate the selection biases affecting 
line-ration selection and demonstrate that the intrinsic 
AGN occupation fraction differs from the the simplest in¬ 
terpretation of the observations, especially for low-mass 
galaxies. 

As discussed in Section 2 and shown in Figure 
[T 1 the well-measured galaxy samples (with line ratio 
(S'/7V)ratio > 3/-\/2) include fewer massive, red, concen¬ 
trated, and low-sSFR galaxies than the parent popula¬ 
tion. In the next subsection we account for AGNs in 
poorly-measured galaxies, showing that the green-valley 
(moderate-sSFR, concentrated) AGN preference is even 
more apparent among the full parent sample of galaxies. 

3.3. The Line-Ratio AGN Fraction Among All Galaxies 

It is only possible to select line-ratio AGNs in galaxies 
with well-measured emission-line ratios. However Figure 
[ 1 ] demonstrates that the well-measured galaxy sample re¬ 
sults in a biased set of galaxy properties. Measuring an 
unbiased AGN fraction requires understanding the likely 
AGN occupation among galaxies with poorly-measured 
{S/N < 3/-\/2) emission-line ratios. 

In Section 3.1 we demonstrated that the V087 AGN 
classification (Equations 2 and 5) tends to select AGNs 
with L[Oiii]/M* > 29.5 (units of erg s“^ and Mq), and 
we impose the same limit to select non-LINER BPT 
AGNs (Equation 6 ). Figure [S] presents the L[Oiii]/M* 
distributions for galaxies with both poorly-measured 
(black) and well-detected emission-line ratios (gray), 
alongside the AGN fraction (red) measured from the 
well-detected galaxy sample. Some of the [Oiii] mea¬ 
surements for poorly-measured galaxies are likely to have 
large error bars, especially at the lowest L[Oiii]/M* ra¬ 
tios. Still, it is clear that most galaxies with poorly- 
measured line ratios typically have log(L[Oiii]/M*) << 
29.5, suggesting that they are unlikely to host AGNs. 
There is only a small tail of poorly-measured galaxies 
with log(L[Oiii]/M*) > 29.5 which might have a signifi¬ 
cant number of “missed” AGNs. 

We account for the potentially missed AGNs by assum¬ 
ing that poorly-measured galaxies have the same AGN 
fraction as a function of L[Oiii]/M, as observed in the 
well-measured galaxies, in each bin of galaxy proper¬ 
ties. Since poorly-measured galaxies have lower typical 
L[Oiii]/M* than well-measured galaxies, bins with sig¬ 
nificant numbers of missed galaxies (i.e., high-mass, red, 
high-sSFR, and concentrated galaxies) are likely to have 
lower AGN fractions than those presented in Section 3.2 
and Figure m for the well-measured galaxy sample. 

Figure | 6 ] presents the total BPT and V087 AGN frac¬ 
tions for the entire parent sample of 317,192 galaxies. 
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Observed AGN Fractions (Well-Measured Galaxies) 

BPT AGN V087 AGN 



Fig. 4.— The observed fractions of BPT AGNs and V087 AGNs among galaxies in the well-measured samples (with (S'/A^)ratio > 3/\/2)? 
in bins of galaxy color, specific star formation rate, and concentration versus stellar mass. Golor is k-corrected to be observed at z = 0.05 
and color—mass bins are tilted (with a slope of 0.5) to match the slope of the blue cloud. Bins with less than 50 galaxies are left blank 
(white). Gray contours represent the overall population of well-measured galaxies. The observed (non-LINER) AGN fraction is highest 
in massive green, moderate-sSFR, and concentrated galaxies and lower on both the low-mass and red (low-sSFR) extremes of the galaxy 
population. 
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BPT Lines V087 + BPT Lines 



Fig. 5.— The distributions of dust-corrected L[Olll]/AF* for 
weak-lined (S/Afratio < 3/\/2) galaxies compared to the AGN 
fraction as a function of L[Olll]/M* measured from strong- 
lined (S’/Aj-atio > 3/\/2) galaxies. The dashed line shows the 
L[Olll]/A^* > 29.5 empirical limit for AGN selection. Most poorly- 
measured galaxies have L[Olll]/AF* < 29.5 and so are typically less 
likely to host AGNs than well-measured galaxies. 

In each bin, we begin with the measured AGN frac¬ 
tion among well-measured galaxies, then add a correc¬ 
tion by assuming the same AGN fraction as a func¬ 
tion of L[Oiil]/M* in the poorly-measured galaxies. 
These corrections are only significant (with >25% as 
many “missed” AGNs as well-measured AGNs) in high- 
mass red, low-sSFR, and concentrated galaxies, and 
the correction in these galaxies results in a smaller 


AGN fraction than found in the well-measured 
pie. The lack of non-LINER AGNs among 

sive quenched galaxies agrees with previous AGb 
studies for line-ratio AGNs (IKauffmann et al.l 

sam- 

mas- 

host 

2003; 

Heckman et al.l 12004 IKauffmann & HeckmanI 

2009; 

Schawinski et al.l 20101: iTanakal 20121), broad-line i 

\GNs 
d X- 

(iTrump et al.l 2013a: IMatsuoka et al.l 1201411. an 

rav AGNs (IMullaney et alJ 20121 Harrison et alJ 120121 

Rosario et al.l l2013allbl: iGhen et al.l 120131: lAzadi et al. 


2015I1 . Otherwise the estimated total BPT and V087 


AGN fractions are very similar to the well-measured 
AGN fractions (Figure S]), and the correction for poorly- 
detected galaxies leads to a qualitatively identical pref¬ 
erence for AGNs to be observed in massive green-valley 
(moderate-sSFR, bulge-dominated) galaxies. 

4. BIASES OF LINE-RATIO GLASSIFIGATION 

The AGN line-ratio signature can be hidden by Hii re¬ 
gion line ratios in galaxies with a significant level of star 
formation. Below we investigate how this “star formation 
dilution” changes with galaxy properties and influences 
the observed AGN fractions. We begin by computing ba¬ 
sic AGN detection limits, quantified by L[Oiii]^qj^/M*, 
for galaxies to host BPT and V087 AGNs. These de¬ 
tection limits implicitly account for the changing spec¬ 
troscopic aperture by using observed Hii region emis¬ 
sion within the same aperture, although we also directly 
investigate aperture effects by comparing the observed 
AGN fractions in different redshift bins. Finally, we 
compute the intrinsic AGN fractions across host galaxy 
properties using Monte Garlo simulations of physically- 
motivated distributions of Eddington ratio and black 
hole / galaxy mass ratios. 

4.1. Basic Limits for AGN Detection 


The biases of BPT and V087 AGN selection can be 
characterized by estimating the threshold at which AGN 
emission exceeds star formation dilution within a 3" 
hber for galaxies of different properties. We quantify 
this threshold using the intrinsic (dust-free) [Oiil] lumi¬ 
nosity of the AGN divided by the galaxy stellar mass: 
L[Oiii]^qj^/M*. This quantity is related to the Edding¬ 
ton ratio Tboi/TvEdd by a [Oiii] bolometric correction 
and the mass ratio For now, we simply as¬ 

sume that L[Oiii]^qi^/M* is monotonically related to 
Fboi/AEdd, reserving a detailed discussion of the func¬ 
tional form until Section 5. 

We estimate the changes in star formation dilution 
with galaxy properties using the typical Hii region line 
strengths in the bins of the various panels in Figure 
m That is, for each bin in color-mass, sSFR-mass, 
and concentration-mass, we calculate the median dust- 
corrected (assuming an intrinsic Ha/H/3 = 2.86) lumi¬ 
nosity of each emission line for galaxies with line ratios 
classified as star-forming in the BPT or V087 diagrams. 
The typical Hii region emission depends both on spectro¬ 
scopic aperture and galaxy properties, and we implicitly 
account for both by using the observed line fluxes of SF 
galaxies in each bin. The AGN detection limit is given 
by the AGN contribution needed to p ush these line ra¬ 
tios on to the AGN/SF diagnostic lines: iKauffmann et al] 
(|2003f) on the BPT (Equation 1), or our empirical line 
on the V087 (Equation 2). For the moment we ne¬ 
glect the AGN/LINER divisions (Equations 5 and 6 ), 
although we include these requirements in our detailed 
models of AGN selection biases in Section 5. We as¬ 
sume a “pure AGN NLR” has log([Oiii]/H/3) = 0.5, 
log([Nii]/HQ:) = 0.0, and Ha/H/3 = 3.1. Because we are 
starting from dust-corrected Hii region emission lines, 
the derived AGN detection limits are similarly dust- 
free, assuming that both NLR and Hii region gas have 
the same dust extinction. We demonstrate that this 
assumption is valid in the Appendix. Using dust-free 
L[Oiii]^qj,j/M* makes it easier to relate to an Edding¬ 
ton ratio, although it does mean that our detection limits 
are not strictly observed quantities. 

Our method of computing AGN detection limits is il¬ 
lustrated by Figure [71 which shows the line ratios and 
strengths (in Lbne/AT*) for two bins in galaxy color- 
mass. Low-mass blue galaxies have higher Lbne/Af* 
from integrated Hii regions than high-mass red galax¬ 
ies. Thus low-mass blue galaxies (which also tend to 
have high-sSFR and low-concentration) require an AGN 
with higher L[Oiii]^Qf.,f/M* to be classified as a BPT or 
V087 AGN compared to massive red (low-sSFR, high- 
concentration) galaxies. 

The L[Oiii]^qj,j/M* detection limits for an AGN to 
exceed the star formation dilution are shown in bins 
of color-mass, sSFR-mass, and concentration-mass in 
Figure jS] At a basic level the detection limits mirror 
the observed AGN fractions: the bins of lowest AGN 
fraction (low-mass, blue, high-sSFR, low-concentration) 
tend to also have the highest AGN det ection thresh- 
olds. If we assume 7 = Mbh/AT* ^ 0.001 ([Haring fc Rixl 
I 200 I assuming M, ^ 1.5Mbuige) and A:[oiii] 100 
([Lamastra et al.ll2009L for log(L[Oiii]^Qf.^) ^ 40), then 
the L[Oiii]^Qf.,[/M* limits translate to log(Lboi/FEdd) 

—2 for blue low-mass (high-sSFR, mixed-concentration) 
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Line Ratio AGN Fractions (All Galaxies) 
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V087 AGN 
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Fig. 6.— Total fractions of observed BPT and V087 AGNs estimated for the entire parent sample of 317,192 galaxies, in bins of 
galaxy color (k-corrected u — z), specific star formation rate, and concentration (-R 90 ,r/-R 50 ,r) versus stellar mass. The AGN fractions are 
determined from the well-measured sample with a correction for the AGNs missed among galaxies with poorly-detected lines, assuming 
the same AGN fraction with L[Oiii]/M* as outlined in Section 3.3. Gray contours show the overall galaxy population from the full parent 
sample. Bins where the correction is >25% of the total AGN fraction are marked with crossed hashes: this only occurs among massive, 
red, low-sSFR, concentrated galaxies. The corrected AGN fractions are very similar to the well-measured AGN fractions in Figure[4| The 
observed BPT and V087 AGN fractions are highest among massive green-valley hosts, with lower AGN fractions for both massive red 
(low-sSFR, concentrated) and low-mass blue (high-sSFR, low-concentration) galaxies. 
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Fig. 7.— An illustration of star formation dilution and AGN detectability for galaxies with different properties. The contours show 
the line ratios and strengths of BPT-classified star-forming galaxies: in blue, low-mass blue galaxies with 8.5 < log(M*/MO) < 8.75 
and —3.275 < {u — z)' — 0.5 log(M*/MO) <= —2.975; in red, high-mass red galaxies with 11 < log(M+/M0) < 11.25 and —2.075 < 
{u — z)' — 0.5log(M*/M©) <= —1.775; and in gray, all galaxies with (S'/A’)ratio > 2>ly/2. Crosses represent the median values of 
the star-forming galaxie s, and filled circles indicate the combined ratios after adding enough “pure AGN NLR” emission to reach the 
IKauffmann et ^ 1|20Q3I) AGN/SF division line. There is more star formation dilution in low-mass blue (high-sSFR, low-concentration) 
galaxies, requiring a relatively more powerful AGN to be selected as a BPT AGN. 


galaxies and log(Lboi/TEdd) ^ —5 for red high-mass 
(low-sSFR, high-concentration) galaxies. Since the 
AGN E ddington ratio distribution i s steeply declin¬ 
ing (e.g., IKauffmann h Heckmanli2009HAird et al .1120 121 : 
iHickox et 'all T^lT ' this behavior leads to a significantly 
lower observed AGN fraction in low-mass galaxies. We 
use more detailed models of Eddington ratio, bolometric 
correction, and MBu/Af* to investigate the AGN selec¬ 
tion biases in Section 5 below. 

Interpreting the biases of BPT and V087 AGN se¬ 
lection crucially d epends on the metric u sed to esti- 
mate AGN power. IKauffmann et ^ (l2003f l (and later, 
iSchawinski et al.l 1201011 nsed L[Oiii] as an indicator of 
AGN strength, concluding that “powerful” (L[Oiii] > 
lO^L©) AGNs are present only in high-mass galaxies. 
These luminous AGNs make up the bulk of the total 
cosmic accretion history of AGNs, since mass accretion 
rate Mbh is related to L[Oiii] by a bolometric correction 
fc[Oiii] and an efficiency 77, Mbh = k[oiii]L[Oiil]/{j]c^)). 
However, the Eddington ratio Tboi/.^Edd is a more effec¬ 
tive indicator of AGN accretion structure and wind/jet 
outflow properties (e.g..lNarav a n fc McClintocldl20n8t iHol 
l2008t iTrump et al.l l2011at iHeckman fc Best! 1201411 . In 
this sense, our Eddington ratio proxy L[Oiil]/M* is an 
effective way to quantify how line-ration selection biases 
relate to AGN-galaxy coevolution. 

The “pure AGN NLR” line ratios assumed here are 
estimated from the locus of AGN-dominated galaxies on 
the BPT diagram. Eigure H] demonstrates that AGN- 
dominated line ratios are most often seen in massive 
hosts. Massive galaxies also tend to have high (near¬ 
solar) metallicities, which means that the pure AGN 
NLR line ratios are appropriate only for metal-rich NLR 
emission. In the Appendix we demonstrate that the ob¬ 
served line-ratio AGN occupation fractions are better de¬ 
scribed by the presence of low-metallicity AGN NLRs in 
low-metallicity galaxies. This effect further worsens the 
AGN detection bias in low-mass galaxies beyond that 


already present in Eigure El 

4.2. Aperture Effects on AGN Detection 

The size of the spectroscopic aperture affects the 
amount of star formation dilution in the observed 
emission-line ratios. Smaller apertnres include less star 
formation re l ative to the nuclear AGN: for example, 
iMoran et al.l (j2002h demonstrated that many nearby 
narrow-line AGNs are hidden by star formation dilution 
in larger spectroscopic apertures. Our objects were all 
observed in 3" fibers with the SDSS spectrograph. Eig¬ 
ure M demonstrates that the 3" aperture covers a dif¬ 
ferent fraction of the galaxy light at different redshifts 
and different galaxy sizes. The median ratio of fiber ra¬ 
dius to r-band half-light radius, Raher/R 50 ,r, ranges from 
^ 0.3 to ^ 0.6 (and Icr bounds extending to ~ 0.2 and 
^ 0.8). The fractional aperture changes are smaller than 
the changes in absolute aperture size because low-mass 
galaxies, which are smaller in physical size, are only ob¬ 
served at low redshift, corresponding to smaller physical 
apertnre sizes. Larger high-mass galaxies are mostly in 
the higher-redshift volume and consequently have larger 
physical aperture sizes. These r-band half-light radii are 
not necessarily equivalent to the emission-line half-light 
radii, but Eigure [9] remains a useful demonstration of the 
changes in aperture effects with redshift. 

The lack of aperture effects with redshift is similarly 
seen in Figure uni which compares observed fractions 
of galaxies classified as BPT AGN and V087 AGN at 
0.01 < z < 0.03, 0.04 < z < 0.06, and 0.08 < z < 
0.1. The “break” from high to low AGN fraction is at 
log(M*) ^ 10 in all three redshift bins. As seen in Figure 
m there are no significant aperture effects causing differ¬ 
ential star formation dilution and AGN identification as 
a function of stellar mass. 

It should be emphasized again that both the AGN lim¬ 
its from Section 4.1 and the Monte Garlo simulations 
of Section 5 implicitly account for aperture effects. We 
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BPT AGN Detection Limits 


V087 AGN Detection Limits 
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Fig. 8.— The detection limits for AGN line ratios to exceed star formation dilution, as a function of galaxy color, specific star formation 
rate, concentration, and stellar mass. Gray contours show the full parent sample of galaxies. The AGN detection limit is quantified as 
L[Oiii]^q]^/M* (in units of erg s“^ Mq^), the dust-corrected AGN [Olll] luminosity divided by the stellar mass. These limits do not 
consider the AGN/LINER classification. The Hii region emission in low-mass blue, high-sSFR, and low-R 9 o,r/R 50 ,r galaxies tends to 
dilute AGN-like line ratios, leading to a L[Oiii]^qj,j/M* detection threshold a factor of ~500 higher than in high-mass red galaxies. 
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Fig. 9.— The fractional aperture covered by the SDSS fiber, 
quantified as the ratio between fiber radius and r-band half-light 
radius, with stellar mass. The thick black line shows the median 
aperture coverage for all galaxies, and the gray dotted lines repre¬ 
sent the 1(T range at each stellar mass. The dashed lines display 
the median fractional aperture in redshift bins of width Az = 0.01 
(i.e., 9 bins from 0.01 < z < 0.02 to 0.09 < z < 0.10). Although 
there are significant aperture changes with redshift, the average 
fractional aperture for the whole sample changes by less than a 
factor of two with stellar mass. 

empirically quantify star formation dilution using line 
fluxes of Hii regions only within the spectroscopic aper¬ 
ture, rather than total-galaxy measurements. Thus even 
if there were aperture effects missed by Figures |9] and 
[TUI they would be accounted for in our models of AGN 
limits and in the resultant intrinsic AGN fraction with 
host properties. 

5. MODELING THE INTRINSIC AGN POPULATION 

In Section 4 we demonstrated that line-ratio AGN se¬ 
lection is biased against host galaxies with low stellar 
mass, blue color, and high sSFR. Such galaxies have more 
“star formation dilution” from optical emission lines pro¬ 
duced in Hii regions. But it remains unclear if the ab¬ 
sence of observed AGNs in low-mass blue disk hosts is 
caused purely by bias, or there is also some change in 
AGN occupation with galaxy properties. 

In this Section we directly model the intrinsic AGN 
population as a function of galaxy properties. Separate 
simulations are run for both the BPT and V087 AGN 
selection methods, comparing the modeled and observed 
fraction of galaxies meeting each set of AGN selection cri¬ 
teria. Within each bin of galaxy properties (color-mass, 
sSFR-mass, or concentration-mass), the simulations fol¬ 
low these basic steps: 

1. A set of “non-AGN” line ratios is randomly drawn 
(with replacement) from the dust-corrected lines 
of observed galaxies with well-measured line ratios 
not classified as BPT AGNs or V087 AGNs. We 
require that these galaxies be face-on (b/a > 0.5) 
and remove outliers in velocity dispersion (using 
iterative 2cr clipping). The number of simulated 
galaxies is the same as the total number of galaxies 
within each bin. 

2. Each simulated “non-AGN” galaxy is assigned an 
AGN with an Eddington ratio Tboi/TEdd randomly 


drawn from a probability distribution. The average 
Eddington ratio is either fixed or allowed to vary 
as a free parameter. 

3. The Eddington ratio is translated to AGN 
emission-line luminosities in three steps. First, 
black hole mass is estimated from galaxy veloc¬ 
ity dispersion (or allowed to be a free parame¬ 
ter). [Oiil] luminosity is then derived from Lboi 
using a (luminosity-dependent) bolometric correc¬ 
tion. From L[Oiii]^Qf,,j, we then compute the other 
line luminosities assuming “pure AGN NLR” line 
ratios. 

4. We compute the number of simulated galaxies with 
resultant AGN+galaxy line luminosities classified 
as BPT AGNs (Equations 1 and 6) or V087 AGNs 
(Equations 2 and 5). The simulated AGN fraction 
is given by the median from 200 Monte Garlo simu¬ 
lations of steps 1-3 above, with an associated error 
given by the normalized median absolute deviation 
(NMAD). 

The assumptions inherent in these steps are motivated 
and discussed in Section 5.1, with detailed tests of the 
assumptions also given in th e Appendix. Inspired by 
the work of lAird et al.l (|2012f) for X-ray AGN hosts, in 
Section 5.2 we test if the observed AGN fractions are 
consistent with uniform SMBH growth regardless of host 
properties. In Sections 5.3 and 5.4, respectively, we in¬ 
stead fit the Eddington ratio and black hole mass as free 
parameters, demonstrating that AGN accretion is corre¬ 
lated with galaxy sSFR and low-mass galaxies may have 
undermassive black holes. 

5.1. Basic Modeling Assumptions 

In each bin of galaxy properties, initial sets of “non- 
AGN” line ratios are empirically drawn from well- 
measured ((<S'/fV)ratio > 3/-\/2) galaxies that have reli¬ 
able velocity dispersions (described below) and do not 
meet the BPT or V087 AGN selection criteria. The 
empirical sets of galaxies implicitly reflect the distribu¬ 
tion of non-AGN line ratios in each bin of galaxy prop¬ 
erties. The breadth of the line-ratio distribution comes 
from the scatter in the mass-metallicity relation (e.g., 
iTremonti et al.l 120041 ) . a range of ionizat ion conditions 
from Hii regions and old star winds (e.g., iSanchez et al.l 
I2015D . and fiber aperture effects for different galaxy sizes 
and redshifts. Beginning the simulations with empirical 
non-AGN line ratios accurately reproduces the observed 
line-ratio distributions (see Section 5.5). Since we do 
not classify LINERs as BPT or V087 AGNs, they are 
included in the “non-AGN” galaxy sets. LINE R galaxies 
are li kely t o host weakly a c cretin g AGNs (e.g., IBB[2(M 
I2009D . but lYan fc Blantonl (I2012D demonstrated that the 
bulk of their emission-line flux is extended rather than 
due to a nuclear AGN. Thus we include LINERs among 
the empirically-drawn “non-AGN” line ratios, simulating 
the addition of more powerful AGNs with Seyfert-like line 
ratios to such galaxies. 

Because the simulations begin with empirical line 
ratios of non-AGN galaxies, they cannot predict the 
AGN fraction in poorly-measured galaxies (which, by 
definition, cannot be reliably classified as star-forming 
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BPT AGN 
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Fig. 10.— The fraction of galaxies classified as BPT AGN and V087 AGN in three bins of redshift, from low to high. Gray contours 
show the distribution of galaxies in each redshift bin. Where galaxies are present in multiple redshift bins, there is little difference in the 
observed AGN fraction, suggesting that AGN selection is minimally affected by the changing aperture. Essentially the increase in galaxy 
size with stellar mass is balanced by a corresponding increase in the typical physical aperture size, since low-mass galaxies are observed 
only at low redshift and high-mass galaxies are mostly at high redshift. 


or AGN). For this reason we only compare (or fit) 
the simulations to the observed AGN fractions among 
well-measured galaxies. The inability to simulate 
poorly-measured line ratios has consequences for mas¬ 
sive red (low-sSFR, concentrated) galaxies, which tend 
to be excluded in the well-measured sample. We in¬ 
vestigate this issue in Section 5.3 by assuming that 
poorly-measured galaxies have the same median ratio of 
Tboi/TEdd/(r[Oiii]jQjg^j/M*) as the well-measured galax¬ 
ies. 

The Eddington ratio distribution of our simulated 
AGNs is parameterized as a Schechter function, following 
iHookins fc HernauistI (|2009ll : 


P(A) 


dt 

dlog A 


= Po 




(7) 


This functional form is effectively a power-law distribu¬ 
tion, with a a “turnover” Eddington ratio of A* = 0.4 
set to impose the Eddington limit A < 1. The Ed¬ 
dington ratio A = Lboi/pEdd, with LEdd — 1-3 x 
10^^{Mbh/Mq) erg s“^, and Pq is a normalization con¬ 
stant. Equation 7 is, at best, an approximation for 
the true Eddington ratio distribution, which is poorly 


known and not well-constrained by observations (e.g., 
iKellv et al.ll201fll l. However, adopting a functional form 
allows us to effectively compare relative Eddington ratios 
across galaxy properties. 

We use a — 0.6 for the power-law slope of 

Equation 7, consistent with the observed Edding¬ 
ton ratio distribution (IHookins fc HernouistJ 120091 
iKauffma nn fc Heck mari|_[2009[l and similar to the dis- 
tribution used bv lAird et al.l (j2012tl . In the Appendix 
we demonstrate that shal lower Eddington ra tio distribu¬ 
tions (e;g;j_^^_R2_Jrom]Hicko^^F^i][20T3 or a = 0.05 
from iSchulze &: Wisotzkill2010fl result in simulations with 
worse fits to the observed AGN fractions. We also pre¬ 
fer a steeper slope in Tboi/TEdd because it qualitatively 
reflects the decreasing effic iency 77 for radiat ively inef¬ 
ficient accretion at low m (|Xie &: Yuanll201^ . With a 
fixed number of Eddington ratios in the random distri¬ 
bution, we set Pq = 1 and numerically control the nor¬ 
malization (and ayerage Eddington ratio) by changing 
the lower bound log(Aniin) of the distribution. 

To conyert Eddington ratio to luminos i ty, we estimate 
black hole mass using the lGiiltekin et all (|2009ll Mbh — cr 
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relation: 

log(MBH/-M^0) = a + /3log(cr/200 km s“^) + eo (8) 


Here a = 8.12, 0 = 4.24, and th e intrinsic scatter 
eo = 0.44. iKormendv &: Hoi (I2013D demonstrate that 
the Mbh — O' relation for classical bulges has signifi¬ 
cantly higher normalization (a = 8.49) and lower s catter 
(eo = 0.29), suggesting that the IGiiltekin et all (120091 1 
normalization is lower due to the inclusion of pseudob¬ 
ulges. However nearly all of our galaxies have insufficient 
data to reliably separat e classical bulg e s and pseudob¬ 
ulges, so we retain the IGiiltekin et al.l (I2009D relation. 
In the Appendix we also test a constant Mbh/M^, ratio, 
finding that it results in a significantly worse fit to the 
data. 

Equation 8 is appropriate only for velocity dispersions 
corresponding to a galaxy bulge. We remove velocity 
dispersions from unresolved rotation or disordered kine¬ 
matics by requiring that the empirically-drawn galax¬ 
ies be face-on {b/a > 0.5), have well-measured veloc¬ 
ity dispersion (error in a <60 km s ~^). and using itera - 
tive 2(7 clipping to remove outliers. iKassin et al.l (l2012f l 
find that most z ^ 0 galaxies with log(M*/M0) > 8.5 
are settled disks, suggesting that unsettled disks with 
dispersion-dominated kinematics are likely to be removed 
by the 2a outlier rejection. The 3” fiber aperture means 
that these velocity dispersions are typically measured 
over a radius smaller than Ve (see Figure |9]). However 
bulges typically have fiat velocity dispers ion profiles over 
0.1 < r/re < 10: iCapoellari et al.l (j2006[) measure a typ¬ 
ical a{r) = tTe(r/re)“° °®® over this range, and so we 
assume a ~ ae for our galaxies. Galaxies with upper 
limits in velocity dispersion are included using survival 
analysis, assuming a normal parent distribution (before 
clipping) described by the median a and its NMAD. 

To convert bolometric luminosity into the extinction- 
corrected [O III] line luminosity we use the following bolo¬ 
metric correction: 


^boi ^ ^^2 ( A[Oiil] 

1040 gj.g g-i yl0"^° erg s“^y 


(9) 


This equation is a power-law fit to the luminosity- 
dependent b o lomet ric corrections found by 
iLamastra et al.l (|2009D . Eq uation 9 is roughly con¬ 
sistent with the findings of iStern fc Laoil (|2012t ) for 
observed (not corrected for extinction) L[Oiil], given the 
1-2 mag of extinction increasing from low to high L[Oiii] 
observed in our galaxies. We adopt an intrinsic scatter 
of 0.4 dex in the L[Oiii] bolometric correction which 
reflects the range of AGN spectral energy distributions 
and details of the emission geometry. This intrinsic 
scatter dominates over t he measurement erro rs in line 
flux and dust correction. IStern fc Laoil (|2012[ 1 measure 
a slightly larger scatter of 0.6 dex for L[Oiii], but some 
of their measured scatter is due to the range of AGN 
extinction. A smaller scatter of 0.4 dex is probably more 
appropriate for the extinction-corrected L[Oiii] used in 
this work. Using a scatter of 0.6 dex in Lboi/A[Oiii] 
instead of the adopted 0.4 dex systematically decreases 
the inferred Eddington ratios by ~0.1 dex, with no 
change in the relative Eddington ratios across galaxy 
properties. 

From all of the (extinction-corrected) 


emission lines in the BPT and V087 diagrams can be 
estimated using line ratios for “pure AGN NLRs”: 

Ha/H/3 = 3.1. (10) 

log([Oiii]/H/3) = 0.5 ± 0.3, (11) 

log([Sil]/Ha) = -0.2 ± 0.2, (12) 

log([Nil]/Ha) = (log([Nii]/Ha)gai) + 0.45 ± 0.2. (13) 

Equation 10 is the intrinsic Balmer decrement 
that best describes typica l AGN NLR conditions 
(jOsterbrock k, Ferlandl[2006[ l. with ionization from the 
AGN causing a value slightly exceeding the Ha/H/3 = 
2.86 of Hii regions. Equations 11 and 12 are empiri¬ 
cally drawn from the AGN locus in the observed BPT 
and V087 diagrams. Equation 13 similarly reflects the 
range of observed AGNs, but also allows for a metallicity- 
dependent AGN NLR. The scatter in each equation was 
chosen to empirically reproduce the observed line-ratio 
distributions in the BPT and V087 diagrams, and effec¬ 
tively accounts for a range in NLR geometry and ioniza¬ 
tion conditions (see Section 5.5). 

To model the metall icity-dependent NLR line ratios, 
we follow the models of iKewlev et al.l (120011 1 , which show 
that lower metallicity decreases the [Nil]/Ha ratio of the 
AGN NLR and galaxy by the same degree, while the 
AGN [Oiil]/H/3 and [Sil]/Ha ratios change little. Simi¬ 
lar metallicity effects on the AGN narrow-line region and 
host galaxy are plausible due to the large (f ew-kpc) size of 
the typical AGN narrow-line region fe.s:.. iBennert et ^ 
l200fR . Most galaxies have fairly fiat metallicity gra- 
dients with a p/H = —0.1 dex/i?e (|Zaritskv et ahUTo^ 

I Sanchez et al.l 1201^ . further suggesting that the NLR 
metallicity should be similar to the galaxy metallicity 
measured within the SDSS fiber (which typically cov¬ 
ers 0.3-0.6i?e, see Figure El- The scatter in Equation 
13 effectively includes the potential effects of scatter in 
metallicity between the galaxy and the AGN NLR. In 
the Appendix we also test a simulation with a uniform- 
metallicity NLR (i.e., a constant NLR [Nil]/Ha ratio), 
finding that it results in a significantly worse fit to the 
observations. Section 5.5 demonstrates that the line ra¬ 
tios and scatter assumed in Equations 10-13 result in 
simulated BPT and V087 diagrams which are similar to 
the observed diagrams. 

5.2. Testing a Uniform Eddington Ratio Distribution 

We begin by testing a “uniform-Aboi /Lk m” model sim¬ 
ilar to that advocated bv lAird et al.l (I2012D : the same Ed¬ 
dington ratio distribution in every bin of galaxy proper¬ 
ties (color, sSFR, and concentration versus stellar mass). 
The simulation uses the Eddington ratio distribution, 
Mbh — cr relation, Lboi/A[Oiii] bolometric correction, 
and metallicity-dependent AGN NLR line ratios outlined 
in Equations 7-13 above. The normalization of Equa¬ 
tion 7 is set to be the same over all galaxy properties: 
log(Aniin) = —5.0 for BPT AGNs and log(Ai„in) = —5.5 

The large size and potential metallicity variation of the 
AGN NLR is in contrast to the much smaller broad-line region, 
which is observed to have s uper-solar metallici t y over a wide range 
in mass and redshift (e.g.. [Nagao et ^120061 : I Juarez et alJr20Q3 : 
IMatsuoka et ani2Qlll'l . 
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BPT AGNs: Uniform L^oi/L^dd. metal-dependent NLR 

Observed BPT AGN Simulated BPT AGN Residuals 
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Stellor Moss [log(M0)] Stellor Moss [log(M0)] Stellor Moss [log(M0)] 


Fig. 11.— The fraction of BPT AGNs with host galaxy color and stellar mass from both observations and the “uniform-Lboi/T/Edd” 
simulation. The simulation uses the assumptions outlined in Equations 7—13: a Schechter function Tboi/T/Edd? relation, luminosity- 

dependent Lboi/T/[Oiii] bolometric correction, and metallicity-dependent AGN NLR line ratios. AGN fraction is indicated by the bin shading 
and gray contours represent the distribution of well-measured galaxies. Similar to the simple AGN limits of Section 4.1 and Figure [S] star 
formation dilution results in fewer AGNs identified in blue galaxies. However the uniform-Lboi/T/Edd model does a poor job of reproducing 
the data, with too many simulated AGNs in red galaxies and too few in high-mass blue/green galaxies. 


V087 AGNs: Uniform Lboi/L^^jj,, metal-dependent NLR 

Observed V087 AGN Simulated V087 AGN Residuals 


N 

I 
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Stellor Moss [log(M0)] Stellor Moss [log(M0)] Stellor Moss [log(M0)] 



Fig. 12.— The fraction of V087 AGNs with host galaxy color and stellar mass from the “uniform-Lbol/T/Edd” simulation compare d to 
the observations. Bin shading indicates the AGN fraction and gray contours represent the well-measured galaxy sample. As in Figure fTTI 
the simulation overpredicts the number of AGNs in massive red galaxies and underpredicts the AGN fraction in massive blue galaxies. The 
uniform-Lboi/T/Edd simulation also fails to reproduce the AGNs in low-mass hosts identified by the V087 but not the BPT. 


for V087 AGNs, with each value set to minimize the 
summed values from comparing to the observations: 


Nbina 


(/obsAGN /simAGN) 
CT^(/obsAGN) 


(14) 


Here /obsAGN and /simAGN are the observed and simu¬ 
lated AGN fractions among well-measured galaxies, and 
gf/nbaAG N) is the binom ial error in /obsAGN (calculated 
following lCameroi]l2011[l . 

The “uniform-Lboi/AEdd” simulation is compared to 
the observed BPT AGN fractions across the color-mass 
diagram in Figure 1111 A uniform Eddington ratio dis¬ 
tribution does a poor job of reproducing the observa¬ 
tions, with significant residuals across the galaxy popu¬ 


lation. The simulations similarly do a similarly poor job 
of reproducing the observed V087 AGN population, as 
seen in Figure [T7] The simulated and observed BPT 
AGN fractions in the sSFR-mass and concentration- 
mass diagrams are also shown in Figures [T3] and fUl 
respectively. In the Appendix we show that the mis¬ 
match with observations is even worse when assuming 
a constant Mbh/-^ 7* ratio or a metal-rich NLR. In all 
cases the “uniform-Lboi/AEdd” model overpredicts the 
observed AGN fraction in massive red, low-sSFR, con¬ 
centrated galaxies, underpredicts in the bluest, highest- 
sSFR, least-concentrated massive galaxies, and overpre¬ 
dicts again in most low-mass galaxies. 

The overprediction of AGNs in massive red (low-sSFR, 
concentrated) galaxies is especially pronounced if we con- 
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BPT AGNs: Uniform Mg^-a, metal-dependent NLR 

Observed BPT AGN Simulated BPT AGN Residuals 
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Fig. 13. — The observed and simulated BPT AGN fractions with host galaxy specific star formation rate and stellar mass, using 
the “uniform-Lboi/TEdd” model. The simulation underpredicts the number of AGNs in moderate-sSFR and moderate-mass galaxies, 
underpredicts the observed AGN fraction at both high and low sSFR and stellar mass. In Sections 5.3 and 5.4 we show that, rather than 
a uniform Eddington ratio, the observations imply a correlation between Eddington ratio and sSFR in massive galaxies, with some hint of 
undermassive black holes in low-mass hosts. 




BPT AGNs: Uniform L^^oi/l-Eddi Mg^-cr, metal-dependent NLR 

Observed BPT AGN Simulated BPT AGN Residuals 
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Fig. 14.— The observed and simulated BPT AGN fractions with host galaxy specific concentration {Rgo,r/R 50 ,r) and stellar mas s, u sing 
“un iform-Lboi/LEdd” model. The uniform Eddington ratio model has the same zebra-like mismatch to the data as in Figures 1111 and 
1131 underpredicting the observed AGN fraction in moderate-mass/moderate-concentration galaxies and overpredicting at low-mass/low- 
concentration and high-mass/high-concentration. 


sider the difference in AGN fractions between the well- 
measured and full samples of galaxies. In Section 3.3 we 
showed that the apparent AGN fraction for massive red 
(low-sSFR, concentrated) galaxies for the full sample is 
lower than in the well-measured sample. Thus the sim¬ 
ulated AGN fraction in such galaxies is likely to be even 
more overpredicted when compared to the full dataset. 

From Figures [TTHT4] we conclude that, with plausi¬ 
ble model assumptions for a bolometric correction, the 
ATbh — O’ relation, and pure AGN line ratios, the ob¬ 
served fractions of line-ratio AGNs are inconsistent with 
a uniform Eddington ratio distribution. In the next 
subsections we instead allow Eddington ratio and black 
hole mass to be free parameters with galaxy properties. 
Rather than a uniform Eddington ratio distribution, the 
observed line-ratio AGN fractions are best reproduced by 


a model where (a) AGN accretion rate is correlated with 
host sSER in massive (log(M*/M0) > 10) host galaxies 
and (b) low-mass (log(M*/M0) < 10) galaxies may have 
a lower black hole occupation function. 

5.3. Non-Uniform AGN Duty Cycle 

Given the poor fit of the “uniform-Lboi/TEdd” simu¬ 
lations, in this subsection we instead allow the normal¬ 
ization of the Eddington ratio (i.e., the average Edding¬ 
ton ratio) to be a free parameter with galaxy proper¬ 
ties. Varying the normalization of the Eddington ra¬ 
tio distribution effectively changes the AGN duty cycle 
(i.e., the number of AGNs above some accretion rate 
threshold). Functionally, this simulation follows steps 
1-4 described in the beginning of the Section. We retain 
the same bolometric correction (Equation 9), Mbh — o’ 
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relation (Equation 8), and metallicity-dependent NLR 
(Equations 10-13) as in the uniform-Lboi/AEdd sim¬ 
ulation. We also use the Schechter function to de¬ 
scribe the AGN Eddington ratio distribution given in 
Equation 7 (with the same declining power-law slope of 
a = 0.6), numerically controlling its normalization us¬ 
ing the lower bound (Lboi/7^Edd)min- However instead 
of a uniform lower bound, we fit over a grid start¬ 
ing from log(Lboi/iEdd)min = -3 and decreasing by 
A log(Lboi/7vEdd)min = 0.05 Until the simulations pro¬ 
duce an AGN fraction equal to the observations. Bins 
with observed AGN fractions consistent with zero (using 
their l-cr binomial error) result in upper limits for the 
Eddington ratio normalization, since lower Eddington ra¬ 
tios would similarly result in zero simulated AGNs. We 
report the best-fit Eddington ratios using the distribu¬ 
tion’s average: for our Schechter function parameteriza- 
tion, (log(Z/bol/7vEdd)) — 0.91 log(Lbol/ 7/Edd)min 0.25. 
A higher average Eddington ratio can also be interpreted 
as a higher AGN duty cycle (i.e. a higher fraction of 
galaxies hosting an AGN above some accretion rate). 

Figures [13] and [12] present the average Eddington ra¬ 
tios from this “non-uniform Lboi/T^Edd” model in bins 
of color-mass, sSFR-mass, and concentration-mass for 
both BPT AGNs and V087 AGNs. The left and mid¬ 
dle columns of each figure show the average Eddington 
ratios calculated from the well-measured set of galax¬ 
ies. In the right column we make a correction to es¬ 
timate the average Eddington ratio for the full galaxy 
sample, assuming that poorly-measured galaxies have 
the same median ratio of Lboi/TEdd/(A[0iii]^Qj3^j/M*) as 
the well-measured galaxies. Within the narrow ranges 
of M* and a in each bin of galaxy properties, this as¬ 
sumption is essentially the same as assuming the same 
fractional AGN contribution to in both well- 

measured and poorly-measured galaxies. Gorrecting for 
poorly-measured galaxies tends to slightly decrease (by 
up to 0.3 dex) the average Eddington ratio in high-mass 
red, low-sSFR, concentrated galaxies (i.e., those galaxies 
most likely to be weak-lined with poorly-constrained line 
ratios). In all three panels, the observed AGN fractions 
in M* < IO^Mq galaxies are frequently consistent with 
zero, resulting in poorly-constrained upper limits for the 
intrinsic Eddington ratios. 

The relationship between AGN Eddington ratio and 
galaxy properties is best understood separately for high- 
mass (log(M*/M0) > 10) and low-mass (log(M*/M0) < 
10) galaxies. We begin with massive galaxies, examining 
the average AGN Eddington ratio as a function of stellar 
mass and specific star formation rate. 

Figure [T3 presents the average Eddington ratio in 
massive (log(M*/M0) > 10) galaxies as a function of 
stellar mass, for both individual bins of sSFR (colored 
points) and the weighted average of the full massive 
galaxy population (black lines). The best-fit line to 
the combined BPT and V087 models has a slope of 
0.065 ± 0.023: consistent with zero. In other words, 
in massive (log(M*/M0) > 10) galaxies, we find that 
AGN accretion is independent of stellar mass. This re¬ 
sult agrees with the uniform £ hni /Lmh with stellar mass 
proposed bv lAird et al.l (|2012l l for X-ray AGN hosts. 

However, Figure [TH| demonstrates that the aver¬ 
age Eddington ratio in massive galaxies is strongly 


non-uniform with galaxy sSFR, increasing by ~2 dex 
with increasing sSFR. Previous work similarly found 
that rapidly accreting AGNs are much more common 
in massive star- forming hosts than in massive quies¬ 
cent hosts le.g.. i Kauffmann et al.l l2003t iHe ckman e t al.l 
2004 liTauffmann fc Heckmanll2009HTrumD et al.ll2013a : 

Rosario et al.ll2013allbl : lMatsuoka et al.ll20l4IAzadi et al.l 

201511 . The solid line in Figure [18] shows the weighted av¬ 

erage Eddington ratio over all mass bins, at a given sSFR. 
We fit a double power-law to the Eddington ratios in each 
sSFR bin for both the BPT and V087 model results. 
With X = log(sSFR/Gyr“i) and y = log(Lboi/TEdd), 
the best-fit (minimum-y^) relationship is given by: 

1. For X = log(sSFR/Gyr“^) < —1.5, 

y = (0.92 ±0.07)x-(2.29 ±0.03). (15) 

2. For X = log(sSFR/Gyr“^) > —1.5, 

y = (-0.10 ± 0.11)a: - (3.79 ± 0.09). (16) 

Error bars in slope and y-intercept are the Ict errors of 
the best-fit values. In other words, Lboi/TEdd increases 
as a power-law with sSFR until log(sSFR/Gyr“^) = 

— 1.5, at which point log(Lboi/TEdd) becomes flat with 
log(sSFR). The slope of the log-log line (exponent 
of the power-law) is well-constrained to be close to 
unity, in striking agreement with the constant ratio of 
Tboi/TEdd to sSFR in th e ensemble of h i gher-r edshift 
X-ray hosts presented by iMullanev et al.l (|2012ll . The 
turnover and lack of dependence of Lboi/TEdd on sSFR 
at log(sSFR/Gyr“^) > —1.5 might be due to AGN con¬ 
tamination in the broad-band photometry at high Ed¬ 
dington ratios: contribution from a blue AGN continuum 
would effectively overe stimate sSFR (e.g.. Appendix A of 
iBongiorno et al.l [20121 1. If sSFR is accurately measured 
and there is negligible AGN contamination, the flatten¬ 
ing of Tboi/TEdd with sSFR may represent a ceiling in 
AGN accretion for host galaxies on the “star forma¬ 
tion mass sequence,” described by log(sSFR /Gyr~^) ~ 

— 1.1 ± 0.3 at z ^ 0 (iWhitaker et al.ll2012h . Similarly, 
iRosario et al.l (|2013bl l and lAzadi et al.l (f2015h argue that 
AGNs are more common in star-forming hosts, but oth¬ 
erwise there is no enhancement of AGN activity with 
additional star formation. 

Interpreting the relationship between AGN accretion 
and galaxy properties from Figures [15] and [T6| is more 
difficult in low-mass galaxies. At the lowest masses 
(log(M*/M0) < 9) the Eddington ratios are essentially 
unconstrained due to zero observed AGNs and extreme 
star formation dilution. Even when the average Edding¬ 
ton ratios are well-constrained at 9 < log(M*/M0) < 10, 
such low-mass galaxies may have black hole masses that 
are not well-described by the Mbh — cr relation. Ap¬ 
plying the Mbh ~ cr relation might overestimate black 
hole mass either because a is poorly measured or results 
from something other than bulge kinematics, or due to a 
black hole occupation function that is less than 100%. In 
the next section we investigate a non-uniform Mbh ~ rr 
relation. 

5.4. Non-Uniform Black Hole Occupation 

iGreene et al.l (|2010fl demonstrated that low-mass 
(^1O®M0) galaxies with accurate megamaser-disk Mbh 
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Best-Fit <Lboi/LEdd> Implied by Observed BPT AON Fraction 

Well-Meosured Goloxles Well-Meosured Goloxles All Goloxies 
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Stellor Moss [log(Mo)] Stellor Moss [log(MQ)] Stellor Moss [log(M0)] 


Fig. 15. — Left: The average Eddington ratios in bins of galaxy color-mass, sSFR-mass, and concentration-mass from the “non-uniform” 
Lboi/^Edd model fits that reproduce the observed BPT AGN fraction for well-measured galaxies in each bin. Gray contours represent the 
well-measured galaxy sample, and hashed bins denote Eddington ratio upper limits where the observed AGN fraction is consistent with 
zero (given its binomial error). Center: The same average Eddington ratios versus stellar mass in each bin, color-coded by {u — 2 )', sSFR, 
or -R 90 ,r/R 50 ,r • Vertical error bars represent the NMAD of the Monte Carlo simulation results. Upper limits due to an observed AGN 
fraction consistent with zero (the hashed bins in the left panel) are shown as thinner lines. Right: Average Eddington ratios for the full 
galaxy sample, assuming that poorly-measured galaxies have the same median L[Oiii]^Qj^/L[Oiii]^Q^g^j as the well-measured galaxies. In 
massive galaxies, the AGN Eddington ratio tends to be highest in hosts with blue color, high sSFR, and low concentration. The connection 
between AGN accretion and galaxy properties is murkier in low-mass hosts. 
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Best-Fit <Lboi/LEdd> Implied by Observed V087 AGN Fraction 


Well-Meosured Goloxles 


Well-Meosured Goloxles 


All Goloxles 



9 10 11 
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9 10 11 
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Fig. 16.— The average Eddington ratios from the “non-uniform” Tbol/TEdd model fits that match the observed V087 AGN fractions 
in bins of galaxy color-mass, sSFR-mass, and concentration—mass. The well-measured galaxy sample is shown by gray contours. Upper 
limits in Edding ton ratio are marked as hashed bins (left panel) and thinner lines (center and right panels). Just as was the case for BPT 
AGNs in Figure ^5] there is a connection between Eddington ratio and host galaxy blue color, sSFR, and low concentration in high-mass 
(log(M*/M 0 ) > 10) galaxies. In low-mass (log(M*/M 0 ) < 10) galaxies the connection between host properties and AGN accretion is 
much less clear. 
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Lboi/LEdd “ Massive Galaxies 



-3 -2 -1 0 

sSFR [log(Gyr"')] 

Fig. 17.— The average AGN Eddington ratio implied by the 
“non-uniform” ibol/^Edd model fits to the observed BPT (top) 
and V087 (bottom) AGN fractions, as a function of galaxy stellar 
mass. Only massive (log(Mt/M 0 ) > 10) galaxies are shown. Bins 
in sSFR are shown by the colored points (with small horizontal off¬ 
sets for clarity) and the black points represent the weighted mean 
at each stellar mass (along with the error of the weighted mean). 
The best-fit line (dashed gray line) is consistent with zero, indicat¬ 
ing that, in massive galaxies, AGN Eddington ratio is independent 
of stellar mass. 

estimates are undermassive by ~0.5 dex with respect to 
the Mbh — O' relation of massive galaxies. This sort 
of non-uniform Mbh ~ t relation, with undermassive 
black holes in low-mass galaxies, can be explained as 
a consequence of black hole seeding mechanisms. Di¬ 
rect collapse of primordial gas is likely t o form the most 
mass i ve (up to -^lO^Mp) SMBH seeds (iBroinm fc LoebI 
l2003t iBegelman et al.l i2006t [Johnson et al.l 120121) . but 
this process becomes ineff icient in low-mass halos 
(|van Wassenhove et al.ll2?)lfll ). Instead, low-mass galax¬ 
ies likely have black holes initially forme d as lower- 
mass (up to ^2OOM|0) Po p III remnants (|Frver et al.l 
120011 iMadau fc Reesll2001l) . Thus low-mass galaxies, 
with low-mass initial seeds and not many mergers over 
their lifetime, might be expected to have a lower black 
hole occupation function com pared to massive galaxies 
(|Volonteri fc Nataraianll2009f) . 

In addition to the physics of black hole seed forma¬ 
tion, it is also possible that problematic a measurements 
might lead to a non-uniform Mbh ~ cr relation. First, 
the measured velocity dispersions at log(M*/M0) < 10 
are frequently below the SDSS instrumental resolution 
(^70 km s“^). Even when well-measured, the velocity 
dispersion of a low-mass galaxy may describe disordered 
kinematics rather than a settled bulge. If the velocity 
dispersion is overestimated for either of these reasons it 
will lead to an overestimate of black hole mass. However, 
the velocity dispersions of low-mass galaxies already im¬ 


Lboi/LEdd “ sSFR in Massive Galaxies 



Fig. 18.— The connection between AGN Eddington ratio and 
host galaxy sSFR for high-mass {log(M*/MQ) > 10.5) galaxies, as 
implied by the BPT AGN (top) and V087 AGN (bottom) frac¬ 
tions. Individual mass bins are shown by the colored points (offset 
horizontally from one another for clarity) and the black points rep¬ 
resent the weighted mean at each sSFR (along with the error of 
the weighted mean). The gray lines are a double power-law fit to 
the combined BPT and V087 data: y = 0.92x — 2.29 at rr < 1.5 
and y = —3.65 at fc > 1.5, where x = log(sSFR/Gyr“^) and 

y = iog(Tboi/iEdd)- 

ply very low-mass black holes: a velocity dispersion of 
35 km s“^, typical of the survival analysis, results in 
only Mbh — W^Mq. This value already approaches the 
“intermediate-mass” black hole regime of direct-collapse 
SMBH seeds. 

Fitting the AGN occupation with Mbh as a free pa¬ 
rameter avoids problematic a measurements and allows 
us to probe the black hole occupation function. We use 
the same approach as our previous simulations, follow¬ 
ing steps 1-4 outlined in Section 5. The normalization of 
the Eddington ratio distribution is fixed by the galaxy 
sSFR according to Equations 14 and 15 above, with 
(log(Tbol/TEdd)) — 0.91 log(Lbol/TEdd)min + 0.25. We 
assume that low-mass galaxies follow the same accretion- 
sSFR relation as massive galaxies, with any deviations 
caused by non-uniform black hole mass. 

Figure [H] and HD] show the best-fit black hole masses 
with galaxy color-mass, sSFR-mass, and concentration- 
mass, implied by both BPT and V087 AGNs. We quan¬ 
tify this Mbh in two ways. The left panels present 
logfMBu/Ma.) , such that bins colored white match the 
iHdring fc Rixl (120041 ) relation with M* = 1.4Mbuige, 
Mbh/M* = 0.001. The center and right panels display 
the deviation from Mbh — t, measured as the log-ratio 
of the best-fit Mbh and the Mbh(o' expected from the 
Mbh — O' relation. Negative log-ratios and redder bins 
denote undermassive black holes, while higher log-ratios 
and bluer bins denote overmassive black holes. Signif¬ 
icantly overmassive black holes are likely due to errors 
in the assumptions, for example, where the Tboi/TEdd“ 
sSFR relation fails to adequately describe the data. 

The “non-uniform Mbh” model does not provide 
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Best-Fit <MgH> Implied by Observed BPT AGN Fraction 


Mbh-M, Devlotion Deviation Deviotion 



8 9 10 11 12 8 9 10 11 12 8 9 10 11 

Stellor Moss [log(M0)] Stellor Moss [log(MQ)] Stellor Moss [log(M0)] 


Fig. 19. — The best-fit average Mbh implied by the observed BPT AGN fraction in bins of galaxy color-mass, sSFR-mass, and 
conc entration-mas s . Lef t: Average black hole mass quantified as log(MBH /-^*)5 where bins colored white match Mbh/-^* = 0.001 
fthe lHaring &: Rd3 l)2004ll relation with M* = 1.4Mbuige)- Gray contours show the distribution of well-measured galaxies, and hashed bins 
denote upper limits where the observed AGN fraction is consistent with zero. Center: Average black hole mass quantified as the deviation 
from Mbh ~ Upper limits with no or very few detected AGNs (the hashed bins at left) are shown as thinner li nes. Right: Th e same 
deviation from Mbh ~ ^ plotted versus stellar mass in each bin, color-coded by {u — z)' , sSFR, or R 90 ,r/R 50 ,r- The lGiiltekin et al.l II2009I) 

Mbh ~ ^ relation and its 1-cr scatter are shown by the dashed and dotted gray lines. Low-mass (10^- IO^^Mq) galaxies may have black 

hole masses that are 0.5—1.5 dex lower than expected from the Mbh ~ ^ relation and their velocity dispersions. 
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Best-Fit <MgH> Implied by Observed V087 AON Fraction 


Mqh-M. Devlotion Deviotion Deviotion 



8 9 10 11 12 8 9 10 11 12 8 9 10 11 

Stellor Moss [log(M0)] Stellor Moss [log(M0)] Stellor Moss [log(M0)] 


Fig. 20. — The best-fit average Mbh^ quantified as log(MBH/-^^*) (left) or deviation from Mbh ~ ^ (center and right), implied by 
the observed V087 AGN fraction in bins of galaxy color-mass, sSFR-mass, and concentration-mass. As before, hashed bins (left panel) 
and thinner lines (center and right panels) denote upper limits where the observed AGN fraction is consistent with zero. Low-mass 

(10®-10^®Mq) galaxies may have slightly undermassive black holes that are ~0.5—1 dex lower than expected from the Mbh ~ ^ relation 

and their velocity dispersions. 
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Observed BPT Simulated BPT 



Fig. 21.— The observed and simulated line ratio distributions in 
th e BPT diagram. The das hed red line represents the empirical line 
of IKauft'mann et al.l 1120031 that we use for A GN/SF classification , 
and the short-dashed blue line also shows the lKewlev et all 11200 11) 
maximal starburst line. The simulated AGN line ratios are broadly 
consistent with the observed data, with minor differences likely 
arising from details of NLR physical conditions. 


Observed V087 Simulated V087 



log([Sll]X6717 + 6733/Ha) log([Sll]A671 7 + 6733/Ha) 

Fig. 22.— The observed and simulated V087 diagrams. The 
dashed red and dotted green lines show our adopted AGN/SF 
and AGN/LINER classification, respectively: we em pirically create 
the AGN/SF line and the AGN/LINER line is from Kewle^^^^ 
(12002). The blue short-dashed line also shows the l^ewIe^^^tT 
11200II) maximal starburst line. As in the BPT diagram in Figure 
1211 there is broad agreement between the observed and simulated 
V087 AGN line ratio distributions. 

any real constraints on Mbh the least-massive 
(log(M*/M0) < 9) galaxies, since their lack of observed 
AGNs and extreme star formation dilution allow almost 
any Mbh ~ cr normalization. The model fits do, how¬ 
ever, provide interesting results for galaxies in the range 
9 < log(M*/M0) < 10. The BPT and V087 AGN frac¬ 
tions both imply undermassive black holes by 0.5-1.5 dex 
in these low-mass hosts. This result is particularly true 
for the least-concentrated and highest-sSFR galaxies, 
which are also most likely to have poorly-understood ve¬ 
locity dispersions (due to disordered kinematics or below 
the instrumental resolution) and/or lower black hole oc¬ 
cupation. We further discuss how the Mbh —o' deviations 
of Figures [19] and [20] relate to black hole occupation and 
the SMBH seed distribution in Section 6.5. 

5.5. Comparing the Observed and Simulated Line Ratio 
Distributions 

The previous two simulations are fit to match the ob¬ 
served AGN fraction as a function of galaxy properties. 
It is also useful to see if they reproduce the observed line 


ratio distribution: this functions as a consistency check 
on our assumed “pure AGN NLR” line ratios. 

Figure [21] compares the observed and simulated BPT 
diagrams, and Figure [52] compares the observed and sim¬ 
ulated V087 diagrams. In both cases the simulated line 
ratios are drawn from a single realization of the Monte 
Carlo simulations in the non-uniform Lboi/FEdd model fit 
to the galaxy color-mass distribution. The line ratio dis¬ 
tributions are similar for fits to bins in galaxy sSFR-mass 
and concentration-mass, and for the non-uniform Mbh 
model, differing only due to the randomly drawn nature 
of each realization. Because the simulations were fit to 
match the observed AGN fractions, the number of galax¬ 
ies classified as AGNs or SF-dominated is identical for 
both observed and simulated line ratios. The simulated 
SF sequence is also nearly identical to the observed ra¬ 
tios, since the simulations begin from a randomly-drawn 
subset of observed non-AGN line ratios. 

Meanwhile the simulated AGN line ratios are broadly 
similar to the observations. The agreement is especially 
good in the [Nil]/Ha and [Sil]/Ha ratios, suggesting that 
our assumed metallicity-dependent NLR model (Equa¬ 
tions 10-13) is reasonably accurate. There are some 
differences in the shapes of the observed and simulated 
[Oiil]/H/3 distributions. In particular, the observed line 
ratios probably do not follow a log-normal distribution, 
but instead depend on physical NL R conditions such as 
ionization or cloud str ucture fe.g.. IGroves et all 12004 
IRichardson et all 120141) . The NLR physical conditions 
may change with AGN Eddington ratio, suggesting that 
it is inaccurate to assume the same normally-distributed 
pure line ratios for every AGN. Still, our simple “evolving 
NLR” description broadly captures the locus and range 
of observed line ratios, and our assumed log-normal dis¬ 
tributions are reasonable in the absence of a better un¬ 
derstanding of AGN NLR physical conditions. 

6. IMPLICATIONS FOR AGN HOST GALAXY PROPERTIES 

The previous sections showed that both the BPT and 
V087 AGN selection methods suffer from significant bias 
in low-mass, blue, high-sSER, low-concentration galaxies 
due to star formation dilution. After carefully modeling 
this bias, we found that the intrinsic AGN population is 
not well-described by a uniform Eddington ratio distribu¬ 
tion over all galaxy properties. Instead there is likely to 
be a connection between AGN accretion and host galaxy 
sSER in massive galaxies. In low-mass galaxies there is 
also some potential for slightly undermassive black holes 
compared to the Mbh — o' relation. 

Below we discuss the ramifications of the line-ration 
selection bias and the implied AGN occupation function 
with galaxy properties. 

6.1. Limitations of Line-Ratio AGN Selection 

As noted in Section 4.1, the biases of line-ratio AGN 
selection depend on how AGN power is quantified. 
The total mass accretion of black holes (integral of 
M over time) is best understand by AGN luminosity, 
since Lboi = pMc^. Meanwhile the Eddington ratio 
Aboi/FEdd gover ns the structure of the accretion flow and 
feedliack mode (iNaravan k, McC lint ockI [2008] : iHol 120081 : 
iTrump et al.ll2011al: iHeckman fc Best]l2014ll 

Eigure [23] presents the fraction of AGNs identified by 
the BPT and V087 methods as a function of both AGN 

























24 


Trump et al. 



37 38 39 40 41 42 43 

log(L[Olll]^GN) [erg/s] 



Fig. 23. — The fraction of simulated AGNs (from the “non- 
uniform Tbol/T^Edd” simulation) identified by the BPT and V087 
AGN classification methods. Detection fractions are presented for 
three different bins of galaxy properties: low-mass and high-sSFR 
(blue lines), high-mass and high-sSFR (green lines), and low-sSFR 
(which also tend to be high-mass, red lines). In the legend, sSFR is 
in units of Gyr“^ and M« is in Mq. In terms of AGN luminosity, 
AGNs are reliably identified at L[Olll] > ® erg s“^ regardless 

of host galaxy properties. The selection function in Tbol/T^Eddi 
on the other hand, is a strong function of galaxy properties, with 
fewer weakly-accreting AGNs identified in low-mass star-forming 
galaxies. Similar trends are seen for both the BPT and V087 clas¬ 
sification methods, although V087 AGN selection has an overall 
lower efficiency. 

luminosity and Eddington ratio. The detection fraction 
is determined from the “non-uniform Tboi/^Edd” sim¬ 
ulation presented in Section 5.3, and is the number of 
simulated AGNs iden tified by the BPT or V087 AGN se¬ 
lection. As found by iKauffiirann et al.l (120031) . the AGN 
detection fraction is high at L[Oiii] > 10^° ® erg s“^ 
(L[Oiii] > Lq) regardless of galaxy properties. In 
other words, the use of line-ratio AGN selection is not 
significantly biased for estimates of total black hole accre¬ 
tion with galaxy properties. However, galaxy properties 
cause a strong bias for line-ratio AGN identification at a 
given accretion rate, similar to that discussed in Section 
4.1 and Figure |S1 The strong “star formation dilution” 
in low-mass star-forming hosts means that, at fixed Ed¬ 
dington ratio, many fewer AGNs are identified than in 
massive and low-sSFR host galaxies. Estimates of AGN 
feedback and fueling efficiency with galaxy properties are 
thus strongly biased when using line-ratio AGN selection. 
Figure!^ also demonstrates that the V087 AGN selec¬ 
tion method has an overall lower efficiency than BPT 
AGN selection. 

The bias against line-ratio AGN selection is worst in 
galaxies with high star formation rates and low metal- 
licity: i.e., the conditions typical of z > 1 galaxies (e.g.. 


iMadau fc Dickinsonll201'i) . Hii regions associated with 
rapid st ar formation may also have harder ionization con¬ 
ditions (iLiu et al.l 20081: iBrinchmumi^Petti ni fc Gharlotl 
l2008t[Seidel et ah 120141 : IShanlev et al.l2015^ . further de- 
creasing the contrast between AGN and star formation. 
Making matters even worse, the BPT loc us at z > 1 is 
suggestive of lower-metall i city AGN NLRs (iKewlev et al.l 
120131 : lJuneau et al.l 12014 iGoil et all I2015D . Thus it is 
likely that Hii regions dominate the emission lines of 
z > 1 galaxies, biasing against AGN line-ratio detection 
to a grea ter degree than i n our low-redshift SDSS sample 
(see also IGoil et al.l[20T5l) . 

Other AGN selection methods that use emission-line 
ratios will be biased in a qualitatively similar fash¬ 
ion to the star formation dilution affecting BPT and 
V087 selection. The w idely-used “mass-excitation” 
(MEx. lJuneau et'n]l201lD and “color-excitation” (CEx, 
lYan et al.ll2011l) methods, which retain [Oiii]/H/3 but re¬ 
place the [Nil]/Ha or [Sii]/Ha line ratios with stellar 
mass or color, will have similar (though more explicit) 
biases against identifying AGNs in blue and low-mass 
galaxies. The CEx method is likely to be more affected 
by lower-metallicity AGN NLRs than MEx AGN selec¬ 
tion, since lower metallicity makes a galaxy bluer without 
changing its mass. 

Line ratios with higher contrast between AGNs and 
Hii regions are likely to be less affected by star for¬ 
mation dilution. High-ionization lines such as [Nev] or 
Hell and (semi-)permitted lines such as Cm] or Civ are 
likely to result in less biased AGN selection since the 
AGN NLR typically has harder ionization and higher 
density than Hii regions. Among ratios of strong lines, 
it is possible th at the [Neiii] /[Oil] ratio used in the 
“TBT” method ([Trouille et alJ[2MTI) offers better diag¬ 
nostic power for AGN selection since [Neiil] has a critical 
density ^10 times higher than [Oiii]. However there are 
many modes of [Neiil] e mission in galaxies, m ost of which 
are not well understood (iZeimann et al.l[2015D . Spatially- 
resolved line ratios also result in higher contrast between 
extended star formation and a nuclear AGN, potentially 
revealing SMBH accretion even when Hii regions dom¬ 
inate the integrated emi ssion lines (e.g., IWright et al.l 
IMntlTrump et ahllWIbl) . 

6.2. AGN Contribution to Emission Lines 

The main goal of this work has been to characterize the 
population of emission-line AGN, including those hid¬ 
den by star formation dilution. But by definition, the 
presence of significant dilution means that Hii regions 
dominate the observed emission lines over AGN NLR 
emission. Many of the AGNs identified in this work have 
little effect on the emission lines from the point of view of 
characterizing galaxy properties (such as star formation 
rate or gas-phase metallicity). 

We use the “non-uniform Lboi/TEdd” model of Section 
5.3 to measure the typical AGN contribution to the ob¬ 
served Ha and [Oiil] emission lines. Quantifying AGN 
luminosity is simpler than inferring the Eddington ratio 
or black hole occupation function, since AGN luminosity 
is directly related to AGN detection above the Hii region 
emission. The assumed Mbh distribution does not actu¬ 
ally matter for determining the typical AGN luminosity: 
for example, a lower black hole occupation would simply 
require an accordingly higher Eddington ratio normal- 
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AGN Fractional Contribution to Ha 

BPT V087 Avg Ha with Moss 
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Fig. 24. — The average fraction of Ho; line emission coming from AGNs in bins of galaxy color-mass, sSFR—mass, and concentration- 
mass. In the left and center panels, gray contours show the full parent sample of galaxies, and bin shading represents the average AGN 
contribution. The right panels show the weighted-average AGN contribution with stellar mass, sSFR, and concentration, weighting by the 
Monte Garlo error in each bin and excluding upper limits. The AGN fraction of Ho emission is highest (up to 15%) in massive red/quiescent 
galaxies, and is lowest (<1%) in low-mass blue/star-forming galaxies. In most galaxies, AGNs contribute <10% to the total Ho emission. 




























































































































































































26 


Trump et al. 


BPT 


AGN Fractional Contribution to 

V087 


[ 0111 ] 

Avg [Olll] with Moss 



0 

0.0 


- 0.5 

- 1.0 

- 1.5 


a: 

m - 2.0 


- 2.5 

- 3.0 


3.5 


S 

01 


3.0 



Avg [Olll] with sSFR 



- 3 . 0 - 2 . 5 - 2 . 0 - 1 . 5 - 1 . 0 - 0.5 0.0 
sSFR log[Gyr"'] 


Avg [Olll] with Cone 



8 9 10 11 12 8 9 10 11 12 

Stellor Moss [log(M0)] Stellor Moss [log(MQ)] 


Cone log[Gyr"'] 


Fig. 25. — The weighted-average frac tion al AGN contribution to the total [Oiii] emission in galaxies binned by color-mass, sSFR—mass, 
and concentration—mass. As in Figure (241 the averages in the right panels are weighted by the Monte Carlo errors in each bin excluding 
upper limits. AGNs are a significant portion (~20-30%) of the average [Oiii] emission in massive red galaxies with low sSFR and high 
concentration, suggesting that the [Oiii] line is not a good indicator of galaxy properties in such systems. 
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Fig. 26.— The average fractional AGN contribution to Ho (top 
panels) and [Olll] (bottom panels), as a function of total line lu¬ 
minosity (left panels) and limiting line luminosity (right panels). 
We account for the line-flux limit of the SDSS in the right panels, 
with the solid line describing the AGN fraction for the magnitude- 
limited (r < 17.77) SDSS and the dotted line using 1/V7iax weight¬ 
ing to es tima te the AGN fraction for a volume-limited survey. As 
in Figure [2^ the AGN contribution to Ho is small (<10%) in most 
galaxies, and only luminous (L(Ho) > 10^^ erg s“^) galaxies are 
AGN-dominated. In contrast, the AGN contribution to [Olll] is 
significant (10-20%) even in weak-lined galaxies and magnitude- 
limited surveys, and galaxies with Z/[Olll] > 10^*^ erg s“^ are dom¬ 
inated by AGNs in [Olll] emission. 


ization to maintain the same number of detected AGNs. 

The fractional AGN contributions to the Hct and [O ill] 
emission lines are shown in Figures [Mj and 1251 respec¬ 
tively. Both the BPT and V087 AGN selection methods 
result in similar estimates of AGN emission-line fractions. 
The average AGN contribution is highest in massive, low- 
sSFR, concentrated galaxies, but is <15% in all galaxies. 
The low average AGN contribution to most galaxies’ Ha 
emission suggests that, for example, Ha-derived star for¬ 
mation rates are generally unbiased by AGN emission in 
similar low-redshift galaxy samples. Meanwhile AGNs 
contribute a significant fraction (^10-30%) of the over¬ 
all [Oiil] emission in massive, low-sSFR, concentrated 
galaxies. In low-mass, high-sSFR, and low-concentration 
galaxies the AGN fraction of both Ha and [Oiii] emis¬ 
sion is negligible (< 1 %), matching our conclusion that 
such galaxies have the largest amount of star formation 
dilution. 

We also measure the fractional AGN contribution as 
a function of line luminosity, shown in the left panels of 
Figure!^ The right panels of Figurepresent the av¬ 
erage emission-line fraction of AGNs for galaxies above 
a given line luminosity limit, useful for estimating the 
typical AGN contribution in emission-line surveys. Each 
cumulative AGN fraction is limited to the redshift range 
at which the the SDSS flux limit (1 x 10“^® erg s“^ cm“^) 
is above the luminosity limit. The solid line is appropri¬ 
ate for the magnitude-limited (r < 17.77) SDSS, while 
the dotted line shows AGN fractions weighted by the 


inverse m aximum det ection volume, 1/I4nax (calculated 
following iHoggI I1999D . giving the AGN fraction for a 
volume-limited survey. The larger number of low-mass 
galaxies (which tend to be star-forming and metal-poor) 
in volume-limited surveys leads to typically lower AGN 
fractions compared to magnitude-limited surveys. In 
general AGNs are an insignificant contribution to Ha 
in all but the most luminous (L(Ha) > 10'^^ erg s“^) 
galaxies and shallowest surveys. However, AGNs are a 
significant fraction of a galaxy’s [Oiii] emission even in 
low-luminosity galaxies and deep magnitude-limited sur¬ 
veys, and AGNs dominate galaxy emission at line lumi¬ 
nosities of L[Oiii] > 10 ^^ erg s“^. 


6 . 3 . AGN Fueling 


Section 5 demonstrated that the Eddington ratio dis¬ 
tributions implied by both BPT and V087 AGNs is 
not uniform across all galaxy properties. At high mass 
(log(M*/M 0 ) > 10), Section 5.3 and Figure [TH| indicated 
that AGN accretion is connected to specific star forma¬ 
tion rate, and otherwise has no dependenc e on stellar 
mass. Here we agree with lAird et al.l (|2012D . who found 
a uniform Eddington ratio distribution with stellar mass 
for X-ray AGN hosts, along with a moderate enhance¬ 
ment of AGN activity with blue galaxy color. The slope 
of unity we derive for the Lboi/LEdd “ sSFR correlation 
at log(sSFR) < —1.5 Gyr“^ is very similar to that of 
IMullanev et al.l (|2012l) for X-ray AGN hosts. We also find 


that the correlation flattens at log(sSFR) > —1.5 Gyr“^. 
This might be caused by AGN contamination affecting 
the estimated sSFR at the highest AGN a ccretion rates. 
It is al so po ssible that, as sug gested by iR.osario et al.l 
(|2013bD and lAzadi et al.! (j2ni5[ ). AGNs prefer hosts on 
the “star-forming mass sequence” but are otherwise un¬ 
correlated with star-forming excess. 

Our estimates of the Eddington ratio distribution are 
less robust at log(M,/M 0 ) < 10, since these galax¬ 
ies have less reliable velocity dispersion measurements 
used to estimate Mbh- Thus we cannot confidently 
claim any non-uniformity in AGN accretion at low stel¬ 
lar masses. In the Appendix we do, however, find that 
the observed BPT and V087 AGN fractions are much 
better described by the iGiiltekin et al.l (|2009D Mbh — t 
relation than by a constant Mbh/AT* ratio. In a galaxy 
of given sSFR, Eddington ratio seems to be uniform 
no matter the bulge fraction, but galaxies with smaller 
bulges tend to have smaller black holes. T his contra¬ 
dicts the constant Mbh/M* ratio used by lAird et ^ 
(|2012D . although some high-redshift observations sug¬ 
gest an evolv ing MBH/Mh,,ire but a constant Mbh/M^, 
fraction fe.g.. [jahnke et al.ll2009t iSchramm fc SilvermanI 


[MllSun e7t al.ll2015l) . 


6 . 4 . AGN Feedback 

The high energy output of AGNs suggests the po¬ 
tential for powerful feedback to quench host galaxy 
star formation. Such feedback could occur ei¬ 
ther through blowout of the star-forming gas by 
radiativelv-driven winds (iSilk fc R^l 19981 lFabianll2002l : 


[Di Matteo, Springel, fc Hernauistll2005D . or by rad io jets 


mechanically heating the gas ( Croton et al.l[200^ . The 
feedback mode likely depends on Eddington ratio, with 
powerful radiative winds likely only at high accretion 
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rates and jets dominating the outflows of weakly ac¬ 
creting AGNs (iNaravan fc Mc Clint ockI l2008t iHol 120081 : 
iTrump et al.ll2011aLlHeckman fc Bestll2014D I3ut in gen¬ 
eral both feedback modes are likely to increase with AGN 
luminosity. 

We define “feedback timescale” as the ratio of a 
galaxy’s gravitational binding energy to the AGN bolo- 
metric luminosity, rfb = t/gai/TAGN- The inverse of this 
timescale quantifies the efficiency of AGN feedback: 

(17) 

Here L 42 = Tagn/( 10'‘^ erg s“i), i?io = it:5o/(10 kpc), 
and Mio = M*/(1O^°M0). Very roughly, at a fixed Ed¬ 
dington ratio and fixed black hole occupation defined 
by Mbh LaGN (Tbol/TEdd)AfBHO'^ (/t)Af*)^ 

(where ft is the bulge-to-total ratio). Meanwhile, galaxy 
mass and size are correlated as M ^ R°‘, with 0.15 < a < 
0.55 i ncreasing from lo w-mass disks to high-mass ellip¬ 
ticals (|Shen et al.ll2003ll . For a uniform Eddington ratio 
distribution this means that AGN feedback efficiency in¬ 
creases with bulge fraction and stellar mass, Tfb ~ 

Of course, the actual situation is somewhat more com¬ 
plicated, as Sections 5 and 6 showed that the Eddington 
ratio distribution is not uniform. 

We use the “non-uniform Tboi/TEdd” simulation of 
Section 5.3 to estimate the AGN feedback efficiency 
across galaxy properties. Just as in Section 6.2, the de¬ 
tails of the Mbh distribution do not matter, since AGN 
luminosity (with respect to the empirically measured star 
formation dilution) is directly constrained by the ob¬ 
served AGN fraction. Our feedback estimates depend 
only on the assumed Lboi/T[Oiii] bolometric correction 
and “pure AGN” line ratios. 

Figure [27] shows the average AGN feedback efficiency 
= TAGN/Hgai in bins of galaxy color-mass, sSFR- 
mass, and concentration-mass. In each bin the feedback 
efficiency is averaged over the Eddington ratio distribu¬ 
tion (which has a declining power-law slope a = 0.6), 
with a small number of high-Lboi/TEdd AGNs and a large 
population of weakly-accreting AGNs. There are slight 
differences in the feedback inferred from the observed 
BPT and V087 AGN populations, largely at low masses 
and high sSFR where the BPT selection results in upper 
limits due to the heavy bias. By either method there 
is little connection between the average AGN feedback 
efficiency and galaxy properties at 2 ; < 0.1: feedback by 
typical AGNs is likely to operate just as efficiently in 
low-mass and high-mass galaxies. The uniformity of 
at high and low sSFR also implies that feedback from 
line-ratio AGNs is not the dominant mode of star for¬ 
mation quenching in most galaxies of any mass or mor¬ 
phology. This is in agreement with recent simulations 
showing that AGN winds, even when present at high ve¬ 
locity, have little ef fect on the larger gas reservoirs within 
their host galaxies ([Gabor fc Bournan3l2014l : lRoos et al.l 

[2?fn . 

6.5. Black Hole Seeds 

The SMBH occupation function in low-redshift galax¬ 
ies carries an imprint of the initial seeding mecha¬ 
nism in the early universe. In particular, low-mass 


(<150 Mq) SMBH seeds from Pop IB remnants are 
likely to be common in most early halos, while mas¬ 
sive (lO'^-lO® Mq) seeds from direct collapse are rare 
in l ow-mass halos due to their higher spin parame¬ 
ter (iBeg elman et al.l 12006]: jyolonteri fc NataraianI 120091 
Ivan Wassenhove et al I TmT Subsequent accretion and 

mergers of these seeds is thoug ht to eventually form the 
local Mbh — a correlation fe.g.. iVolonteri fc Starkl]2Ml] : 
lVolonteri]l2012ll . Due to the high direct-collapse seed¬ 
ing rate and high merger rate, essentially every massive 
galaxy is likely to end up on the Mbh — cr relation. Both 
the merger rate and the fraction of massive seeds are 
thought to be lower for low-mass (log(M*/M 0 ) < 10 ) 
hosts, and so many of these galaxies may have under- 
massive black holes compared to the Mbh ~ cr relation. 

We determine the “SMBH occupation function” as the 
fraction of galaxies with black hole masses well-described 
by the Mbh ~ cr relation. Our galaxies are assumed to 
contain a bimodal population of black holes, with one 
peak of massive seeds described by the iGiiltekin et al.l 
1200 ^ Mbh — a relation (Equation 8 ), and a second log¬ 
normal distribution of lower-mass Pop IB seeds centered 
at 150Mq ±0.5 dex. From this distribution the SMBH 
occupation function is fit using a simulation similar to 
that in Section 5.4, with the same bolometric correction 
(Equation 9), metallicity-dependent NLR (Equations 10- 
13), and sSFR-dependent Eddington ratio distribution 
(Equations 14 and 15), while using the bimodal distri¬ 
bution for Mbh- From this distribution, the SMBH oc¬ 
cupation function is described by the fraction of galaxies 
hosting black holes drawn from the higher-mass Mbh — o' 
peak. 

Figure [28l shows the SMBH occupation function from 
our model fits to the observed BPT and V087 AGN frac¬ 
tions. There are some differences between the two selec¬ 
tion methods, likely due to the uncertainties associated 
with NLR metallicity in low-mass galaxies. In addition, 
the velocity dispersions in log(M*/M 0 ) < 10 galaxies are 
poorly understood, due to both the low SDSS resolution 
(^70 km/s) and the difficulty in interpretation (bulges 
versus disordered disks). Nonetheless both the BPT and 
V087 AGN fractions imply a slightly lower black hole 
occupation of ^30-50% in low-mass and disk-dominated 
hosts. This res ult is roughly consi stent with the indepen¬ 
dent study of iMiller et al.l ([201511 . who similarly found 
a marginally lower SMBH occupation function of a few 
tens of percent for \og{M^/ Mq) < 10 X-ray AGN hosts. 
Confirming this marginal evidence and robustly deter¬ 
mining the SMBH black hole occupation function at low 
stellar mass requires higher resolution spectroscopy than 
available from our SDSS sample. 

7. SUMMARY 

We constructed a sample of over 300,000 galaxies from 
the Sloan Digital Sky Survey to investigate the biases of 
line-ratio AGN selection and recover the intrinsic AGN 
population across a range of galaxy properties. While 
line-ratio AGNs are observed to be most common in mas¬ 
sive green-valley hosts, we demonstrated that this result 
is a selection effect caused by the bias from “star for¬ 
mation dilution” in low-mass and star-forming galaxies. 
After accounting for this bias, we find that AGNs are 
most common in massive galaxies with high specific star 
formation rates, implying that SMBH accretion and star 
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Fig. 27. — The average AGN feedback efficiency = TAGN/C^gal with galaxy color-mass, sSFR-mass, and concentration-mass, inferred 
by both the BPT and V087 AGN selection. The average feedback efficiency is given by the bin shading in the left and center panels, with 
the full parent sample of galaxies shown by gray contours. The feedback efficiency is shown with stellar mass, sSFR, and concentration 
in the right panels, computed as weighted averages using the Monte Carlo error in each bin and excluding upper limits. AGN feedback 
efficiency is fairly uniform with galaxy properties, with no evidence that AGNs dominate the quenching of star formation for most galaxies 
of any type. 
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Black Hole Occupation Function 
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Fig. 28. — The SMBH occupation function derived from our model fit to the observed BPT AGN and V087 AGN fractions. Here we 
assume that the black hole population is bimodal with two log-normal peaks, one from Pop III star remnants at ISOMq (with ±0.5 dex 
scatter) and one from massive seed “fossils” described by the lGiiltekin et al.l II2009I1 Mbh relation (Equation 14). With the caveat that 
a may be inaccurate (in both measurement and interpretation as bulge kinematics) in our low-mass galaxies, there is some hint of a lower 
SMBH occupation function in log(M,/M0) < 10 hosts. 
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formation are fueled by the same gas reservoir. AGNs 
contribute little to the overall Ha and [Oiii] emission 
lines in low-mass galaxies, but their feedback effects are 
likely to be just as efficient at all stellar masses and host 
morphologies. There is marginal evidence that the black 
hole occupation function may be a factor of a few lower 
in \og{A'U/M q) <10 hosts, although this result is not 
robust due to poorly-understood velocity dispersions in 
these galaxies. Higher resolution spectroscopy is needed 
to better constrain SMBH occupation function, while se¬ 
lection methods via spatially resolved line ratios or high- 
ionization lines suffer less bias and can better reveal the 
connection between AGN accretion and star formation. 
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APPENDIX 

TESTING OTHER MODELING ASSUMPTIONS 

In this Appendix we justify the modeling assumptions used in the simulations of Section 5 (and outlined in Equations 
7-13) by testing other plausible assumptions. The models used the same dust extinction for both Hii and AGN NLR 
emission lines, and we test that assumption here. The original “uniform-Lboi/AEdd” simulation is also compared with 
additional simulations that replace the metallicity-dependent NLR ratios with metal-rich (constant) NLR ratios, or 
replace the Mbh — <7 relation with a constant Mbh/A 7* ratio. We also test different power-law slopes for the Schechter 
function describing the Eddington ratio distribution. 

Dust Corrections 

Our estimates of the intrinsic AGN occupation fraction used dust-free emission lines, beginning from dust-corrected 
emission lines observed for star-forming galaxies, and then adding dust-free model AGN NLR emission lines. Thus the 
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Dust Properties for Goloxies ond BPT AGNs 
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Fig. 29.— The distribution of dust extinction for galaxies classified as BPT AGNs and non-AGNs (SF galaxies and LINERs). We 
quantify the dust extinction as AlogfHa/H/^) = log(HQ;/H/3 )r»hs — log(Ho/H/3)i, where the intrinsic Balmer decrement (Ho/H^)! is 3.1 for 
AGNs and 2.86 for non-AGNs flOsterbrock Ferlandll2006l) . AGNs have marginally more dust extinction than the full set of non-AGNs 
(left panel), but this is because BPT AGNs are identified only in massive galaxies. AGNs have nearly identical dust extinction compared 
to non-AGNs of the same stellar mass (right panel). 
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Fig. 30.— The observed and simulated fractions of BPT AGNs with host galaxy color and stellar mass for the “uniform Tboi/TEdd” 
simulation with metal-rich NLR line ratios. Bin shading indicates the AGN fraction and the gray contours represent the distribution of 
well-measured galaxies. Using a metal-rich NLR results in a larger overprediction of BPT AGNs in low-mass galaxies compared to the 
metallicity-dependent NLR used in Section 5.2. 


TABLE 1 

X^/1000 Values for Simulated AGN Fractions 


Model 

(u — z)' — Mt 

BPT 

sSFR -M, 

Cr-M, 

1 

1 

V087 
sSFR -M* 

Cr - M, 

Uniform Lboi/TEdd? ATbh — cr? metal-dependent NLR 

12.8 

10.7 

9.83 

2.62 

2.32 

2.37 

Uniform Lboi/TEdd^ ^BH — metal-rich NLR 

23.1 

19.2 

15.7 

2.62 

2.32 

2.37 

Uniform Lboi/TEdd? -^BH ~ -^*5 metal-rich NLR 

46.0 

33.2 

34.0 

4.59 

4.34 

4.97 

Shallow (CK = 0.2) Lboi/LEdd distribution 

24.5 

18.6 

16.7 

4.59 

4.34 

4.97 

Flat (o = 0.05) Lboi/^Edd distribution 

28.2 

21.4 

18.9 

5.48 

4.09 

4.27 
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BPT AGNs: Uniform metal-rich NLR 

Observed BPT AGN Simulated BPT AGN Residuals 
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Fig. 31.— The observed and simulated fractions of BPT AGNs with host galaxy color and stellar mass for the “uniform Tboi/T/Edd” 
simulation with metal-rich NLR line ratios and a constant Mbh/^* ratio. As before, the AGN fraction is indicated by the bin shading and 
the well-measured galaxy distribution is given by the gray contours. Assuming a co nsta nt MsH/Af* ratio results in a significantly worse 
match to the observations than the simulation using the Mbh ~ ^ relation in Figure !^ (and in Section 5.2). 
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Fig. 32.— A comparison of the observed and simulated BPT AGN fractions in the galaxy color-mass diagram, using a simulation with 
a “shallow” Eddington ratio distribution of power-law slope a = 0.2. The shallowLuol/T/Edd distribution results in larger residuals than 
the distribution with a = 0.6 adopted in the remainder of this work (e.g., Figure [11} . 

models implicitly assume the same dust extinction for both star-forming galaxies and AGNs: a plausible assumption, 
given that the AGN NLR gas is at >kpc scales from the galaxy center ( Bennert et al.l[2002f) . similar to locations of 
Hii regions. We directly compare the dust extinction of BPT-classified AGNs and non-AGNs in Figure (291 quantified 
as the excess Balmer decrement Alog(Ha/H/3) = log(Ha;/H/3)nbs — loK(Ha/H/3)i. Here the int rinsic Balmer decrement 
(Ha/H/3)i = 2.86 for Hii regions and (Ha/H/3)i = 3.1 for AGNs (|OstCTbrock fc Ferlandll20d^ . Gompared to galaxies 
of the same (high) stellar mass, AGNs have nearly identical Balmer decrements, justifying our assumption that both 
Hii region and AGN NLR gas are affected by the same dust extinction. 

Constant-Metallicity NLR and Constant Mbh/M* Ratio 
To test a constant, metal-rich NLR we replace Equation 13 with: 

log([Nii]/Ha) = 0.0 ± 0.2, (Al) 

In other words, we fix the AGN [Nil]/Ha ratio to a normal distribution rather than having it depend on the [Nil]/Ha 
ratio of the host galaxy. The metal-rich NLR assumption uses the same equations as the Section 5 simulations 
(Equations 10-12) for the other AGN NLR line ratios: this means there is no difference between the metallicity- 
dependent and metal-rich V087 AGN simulations. We follow the same steps outlined in Section 5.2 to create a 
“uniform-Lboi/TEdd” simulation with a metal-rich AGN NLR, using the same log(Ainin) = —5.0 for BPT AGNs and 
log(Aniin) = —5.5 for V087 AGNs to normalize the Eddington ratio distribution and minimize the total in comparing 
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BPT AGNs: Flat a = 0.05 Edd ratio distribution 

Observed BPT AGN Simulated BPT AGN Residuals 
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Fig. 33.— A comparison of the observed and simulated BPT AGN fractions in the galaxy color-mass diagram for a simulation with a 
“fiat” Eddington ratio distribution of power-law slope a = 0.05. The simulation with a flat Tboi/^Edd distribution produces a significantly 
poorer fit to the data than the steeper a. = 0.6 Tbol/-^Edd distribution. 

to observations. This uniform-Lboi/AEdd, metal-rich NLR simulation for BPT AGNs is compared to observations in 
the color-mass diagram in Figure 1301 

We additionally test a different relationship between black hole mass and galaxy properties by replacing Equation 8 
with; 

Mbh/M* = 0.001, (A2) 

with an intrinsic s catter of 0.5 dex. T his is consistent with the A/Bn/ATbuige relation of lHaring fc Rixl (120041) . assuming 
Mbuige = Af*/1.4. lAird et HI (120121) used a similar constant Mbh/AT* ratio in their study showing that X-ray AGNs 
have a uniform Eddington ratio distribution over a wide range of host galaxy stellar mass. In the local universe ve locity 
dispersion correlates with black hole mass much better than does total stellar mass ('e.g.. iKormendv fc Holl2013ll. but 
several z > 1 studies suggest that Mnn/AIbiiiyR ev olves with redshift while Mbh/M^, is constant (e.g.. iJahnke et ahl 
l2009t[^hramm &: Silvermanll2013llSun et al.ll2015|l . Thus there is some motivation for comparing a constant Mbh/AI* 
ratio with the Mbh ~ c assumption used in Section 5. We create a “uniform-Lboi/AEdd” simulation with both a 
metal-rich NLR and a constant Mbh/AI* ratio following the same steps as in Section 5.2, minimizing the total 
by setting log(Aniin) = —5.75 in both the BPT and V087 simulations. Eigure |3T] displays this uniform-Lboi/AEdd, 
metal-rich NLR, constant-MBH/.^if* simulation for BPT AGNs in the color-mass diagram. 

Table [T] compares the y^/lOOO values from each simulation for BPT and V087 AGN fractions in the color-mass, 
sSER-mass, and concentration-mass host galaxy diagrams. As seen in Eigures [30] and [311 using a metal-rich AGN 
NLR or a constant Mbh/A 7* relation results in significantly worse matches between the simulated and observed AGN 
fractions. Thus we are justified in using the metallicity-dependent NLR and Mbh — <2 ' relation for our simulations in 
Section 5. 


Different Power-Law Slopes for the Eddington Ratio Distribution 

The simulations in Sections 5 and 6 used a Schechter form of the Eddington ratio distribution (Equation 7) with 
a power-law slope a = 0.6. Observati ons of X-ray and narrow-line AGNs, as well as models for AGN fueling, ar e 
consistent with this parameterization ([Hopkins k, HernauistI 120091 : iKauffmann fc HeckmanI l2009t lAird et all l2012f) , 
but neither are strongly co nstraining du e to th eir considerable uncertainties. Other observations suggest a somewhat 
shallower power-law slope: iHickox et al.1 (j2014D argue that a = 0.2 better describes the observed LagnSFR correla¬ 
tion, and lSchulze fc Wisotzkil ([201011 measure a = 0.05 for low-redshift broad-line AGNs. Here we test the effects of 
adopting shallower Eddington ratio distributions on the fit between the simulated and observed AGN fractions for the 
uniform-Lboi/AEdd model. 

Eigure [32l compares the uniform-Lboi/AEdd simulation with an Eddington ratio distribution of slope a = 0.2 with 
the observed AGN fraction. This “shallow” Lboi/AEdd distribution was normalized (minimizing the y^) using a lower 
log(Aniin) = —7 for BPT AGNs and log(Aniin) = —6.5 for V087 AGNs. The fit to the observed AGN fractions is 
significantly worse than the fiducial a = 0.6 used in Sections 5 and 6, with y^ values ^1.5x higher, as shown in Table 

m 

The uniform-Lboi/AEdd simulation with a “flat” a = 0.05 Eddington ratio distribution is compared to the obser¬ 
vations in Eigure [33l Normalizing the flat Lboi/AEdd distribution requires a very low minimum Eddington ratio: 
log(Aniin) = —10 for BPT AGNs and log(Ai„in) = —9 for V087 AGNs. The y^ values for the a = 0.05 distribution 
are roughly double those of the fiducial a = 0.6 Eddington ratio distribution (see Table [1]), with a higher-amplitude 
“striping” pattern of residuals in the color-mass diagram. Eor a uniform Eddington ratio distribution in all host galaxy 
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types, the steeper a = 0.6 Eddington ratio distribution provides the best (lowest-y^) match to the observed AGN 
fractions. 



